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Abstract 

The mass and width of the W boson are measured using e+e" — W+W~ events from 
the data sample collected by the OPAL experiment at LEP at centre-of-mass energies 
between 170 GeV and 209 GeV. The mass (mw) and width (Fw) are determined using 
direct reconstruction of the kinematics of W+W~ qqliy and W+W~ — > qqqq events. 
When combined with previous OPAL measurements using W^W^ Iviiy events and 
the dependence on mw of the WW production cross-section at threshold, the results are 
determined to be 

TOW = 80.415 ± 0.042 ±0.030 ±0.009 GeV 
Fw = 1.996 ±0.096 ±0.102 ± 0.003 GeV 

where the first error is statistical, the second systematic and the third due to uncertainties 
in the value of the LEP beam energy. By measuring tow with several different jet algo- 
rithms in the qqqq channel, a limit is also obtained on possible final-state interactions due 
to colour reconnection effects in W+W" — > qqqq events. The consistency of the results 
for the W mass and width with those inferred from other electroweak parameters provides 
an important test of the Standard Model of electroweak interactions. 
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1 Introduction 

The measurement of the mass of the W boson (mw) is one of the principal goals of the physics 
programme undertaken with the LEP e+e~ collider at CERN. Within the Standard Model of 
electroweak interactions, the W mass can be inferred indirectly from precision measurements 
of electroweak observables, in particular from e+e" — > Z events at centre-of-mass energies 
(v^) close to the peak of the Z resonance (around 91 GeV), studied extensively at LEPl and 
SLD [1]. These measurements currently give a prediction for mw with an uncertainty of 
32 McV, or 23 MeV if the measurement of the mass of the top quark from the Tevatron [2] 
is also taken into account. Direct measurements of the W mass with a similar precision are 
therefore of great interest, both to test the consistency of the Standard Model and better to 
constrain its parameters (for example the mass of the so-far unobserved Higgs boson), and 
to look for deviations signalling the possible presence of new physics beyond the Standard 
Model. Such measurements became possible at LEP once the centre-of-mass energy was 
raised above 160 GeV in 1996, allowing the production of pairs of W bosons in the reaction 
e''"e~ — W'*"W~. Measurements of the width of the W boson (Fw) can also be carried out 
at LEP, providing a further test of the consistency of the Standard Model. 

This paper presents the final OPAL measurement of the mass and width of the W boson, 
using direct reconstruction of the two boson masses in e"''e~ W+W~ — > c^qiv and qqqq 
events recorded at e"*'e~ collision energies between 170 GeV and 209 GeV. The result for 
mw is combined with a measurement using direct reconstruction in the Ivlv final state [3] 
and a measurement from the dependence of the WW production cross-section on mw at 
161 GeV [4]. This paper supersedes our previous results [5, 6, 7] obtained from the data 
with =170-189 GeV. 

Three methods are used in this paper to extract mw, all based on similar kinematic fits to 
the reconstructed jets and leptons in each event. The principal method, the convolution fit, 
is based on an event-by-event convolution of a resolution function, describing the consistency 
of the event kinematics with various W boson mass hypotheses, with a Breit-Wigner physics 
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function dependent on the assumed true W boson mass and width. The convolution fit 
is used to obtain the central results of this paper, but is complemented by two other fit 
methods of slightly lower statistical precision: a reweighting fit based on fitting Monte Carlo 
template distributions with varying assumed W mass and width to the reconstructed data 
distributions, and a simple analytic Breit-Wigner fit to the distribution of reconstructed 
W boson masses in the data. Complete analyses, including systematic uncertainties, have 
been performed for all three methods, providing valuable cross-checks of all stages of the 
analysis procedure. The convolution and reweighting fits also measure the W width; the 
convolution fit is again used for the central results, and the reweighting fit provides a cross- 
check including all systematic uncertainties. The Breit-Wigner fit does not measure the W 
width, but an additional independent convolution-based method is used to provide a second 
statistical cross-check in the qq£v channel. 

The dominant systematic error in the qqqq channel comes from possible final-state in- 
teractions (colour rcconnection and Bose-Einstein correlations) between the decay products 
of the two hadronically decaying W bosons. According to present phenomenological models, 
these interactions mainly affect soft particles, and the uncertainties can be reduced by re- 
moving or deweighting soft particles when estimating the directions of jets. Such a method 
is used for the qqqq channel measurements of mw from all three fit methods in this paper. 
Conversely, the effect of final-state interactions can be enhanced by giving increased weight 
to soft particles, and this is used to place constraints on possible colour reconnection effects. 

This paper is organised as follows. The OPAL detector, data and Monte Carlo samples 
arc introduced in Section 2, followed by a brief description of the event selection in Section 3. 
Elements of the event reconstruction and kinematic fitting common to all three analysis meth- 
ods are discussed in Section 4, followed by a detailed description of the individual convolution, 
reweighting and Breit-Wigner fits in Sections 5-7. Systematic uncertainties, which are largely 
common to all three methods, are described in Section 8. Finally the results are summarised 
in Sections 9 and 10. 

2 Data and Monte Carlo samples 

A detailed description of the OPAL detector can be found elsewhere [8]. Tracking of charged 
particles was performed by a central detector, enclosed in a solenoid which provided a uniform 
axial magnetic field of 0.435 T. The central detector consisted of a two-layer silicon microvertex 
detector, a high precision vertex chamber with both axial and stereo wire layers, a large 
volume jet chamber providing both tracking and ionisation energy loss information, and 
additional chambers to measure the z coordinate of tracks as they left the central detector.^ 
Together these detectors provided tracking coverage for polar angles |cos^| < 0.96, with a 
typical transverse momentum (px) resolution^ of CFp^/p^ = i/(0.02)^ -|- (O.OOlSpr)^ with pr 
measured in GeV. The solenoid coil was surrounded by a time-of-flight counter array and a 
barrel lead- glass electromagnetic calorimeter with a prcsamplcr. Including also the endcap 
electromagnetic calorimeters, the lead-glass blocks covered the range jcos0| < 0.98 with a 
granularity of about 2.3° in both 6 and <p. Outside the electromagnetic calorimetry, the 

right handed coordinate system is used, with positive z along the e~ beam direction and x pointing 
towards the centre of the LEP ring. The polar and azimuthal angles are denoted by 6 and 0, and the origin 
is taken to be the centre of the detector. 

^The convention c = 1 is used throughout this paper. 
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Year 


■^s range 




/Ldt 


qqei/ 




qqrv 


qqqq 




(GeV) 


(GcV) 


(pb-i) 


obs. 


exp. 


obs. 


exp. 


obs. 


exp. 


obs. 


exp. 


1996 


170-173 


172.1 


10.4 


22 


20 


15 


19 


13 


10 


60 


58 


1997 


181-184 


182.7 


57.4 


134 


122 


117 


124 


118 


124 


437 


446 


1998 


188-189 


188.6 


183.1 


388 


413 


422 


417 


444 


425 


1551 


1511 


1999 


192-202 


197.4 


218.5 


524 


512 


489 


518 


559 


526 


1924 


1891 


2000 


200-209 


206.0 


219.6 


506 


524 


530 


525 


555 


543 


1921 


1925 


Total 


170-209 


196.2 


688.9 


1574 


1591 


1573 


1603 


1689 


1628 


5893 


5831 


Estimated selection efRciency (%) 




85 




89 




68 




86 


Estimated purity (%) 






92 




92 




73 




79 



Table 1: Observed and expected numbers of candidate WW events, together with the collision 

energy range, mean energy and integrated luminosity, in each year of data taking. The 
efficiencies and purities of the event selections, estimated from Monte Carlo events, are also 
given. 

magnet return yoke was instrumented with streamer tubes to form a hadronic calorimeter, 
with angular coverage in the range | cos^| < 0.91 and a granularity of about 5° in 9 and 7.5° 
in (f). The region 0.91 < | cos 9\ < 0.99 was instrumented with an additional pole-tip hadronic 
calorimeter using multi-wire chambers, having a granularity of about 4° in 9 and 11° in (p. 
The detector was completed with muon detectors outside the magnet return yoke. These were 
composed of drift chambers in the barrel region and limited streamer tubes in the endcaps, 
and together covered 93 % of the full solid angle. The integrated luminosity was evaluated 
using small angle Bhabha scattering events observed in the forward calorimeters [9]. 

The data used for this analysis were taken at centre-of-mass energies between 170 GeV 
and 209 GeV during the LEP2 running period from 1996 to 2000, and correspond to a total 
integrated luminosity of about 689 pb~^. In the year 2000, LEP was operated in a mode 
where the beam energy was increased in ~ 0.5 GeV steps during data taking several times 
in each collider fill. Data taken during these 'miniramps' (approximately 1 % of the total 
year 2000 data sample) are excluded from the analysis as the beam energy is not precisely 
known. A detailed breakdown of the energy ranges and integrated luminosities in each year 
of data taking is given in Table 1. In addition, e"'"e~ — > Z events recorded at ^/s « 91 GeV 
were used to calibrate the leptonic and hadronic energy scales and to study the modelling 
of the detector response by the Monte Carlo simulation. These events were recorded during 
dedicated runs at the beginning of each year, and also at intervals later in the data-taking 
periods to monitor the stability of the detector performance with time. They amount to 
a total integrated luminosity of about 13pb~^, corresponding for example to about 400 000 
hadronic Z decays. 

Large samples of Monte Carlo simulated events have been generated to optimise and 
calibrate the W mass and width analysis methods, and to study systematic uncertainties. 
The relevant contributions to the e+e~ — »• qq£u and qqqq topologies studied in this paper can 
be divided into four-fermion and two-fermion processes [10]. As defined here, four-fermion 
final states (e^e~ 4f) include contributions from both e^e~ W'''W~ 4f and e'''e~ — >■ 
ZZ 4f, but exclude multi-peripheral diagrams resulting from two-photon interactions, 
which have a negligible probability of being selected by the analysis requirements and are 
not considered further. Most four-fermion final states were simulated using the KoralW 
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1.42 program [11], which uses matrix elements calculated with grc4f 2.0 [12]. These samples 
were split into two parts, corresponding to four-fermion final states which could have been 
produced from diagrams involving at least one W boson (referred to collectively as WW events 
below), and others (referred to as ZZ events, but including some diagrams not involving two 
Z bosons). Most WW events were generated with mw = 80.33 GeV and Fw = 2.09 GeV, 
but samples with other W masses and widths were also produced in order to calibrate and 
test the fitting procedures. The running width scheme for the Breit-Wigner distribution as 
implemented in KoralW was used throughout. Four-fermion background from the process 
e+e~ — >■ e+e^qq (included in the ZZ sample) was simulated using grc4f. The only important 
two-fermion background process is e+e" — > Z/7 — > qq, generated using KK2f 4.13 [13], with 
Pythia 6.125 [14] as an alternative. 

Hadronisation of final states involving quarks was performed using the Jetset 7.4 model 
[15], with parameters tuned by OPAL to describe global event shape and particle production 
data at the Z resonance [16]. This hadronisation model and parameter set is denoted by JT. 
To study systematic uncertainties related to hadronisation, the same two- and four-fermion 
events have been hadronised with various alternative hadronisation models and parameter 
sets: Jetset 7.4 with an earlier OPAL-tuned parameter set based primarily on event shapes 
[17] (denoted JT'), Ariadne 4.08 [18] with parameters tuned to ALEPH data [19] (denoted 
by AR), Ariadne 4.11 (AR') and Herwig 6.2 [20] (HW), both with parameters tuned to 
OPAL data. The possible effects of final-state interactions in e+e" — >■ W+W" qqqq events 
have been studied using colour reconnection models implemented in Pythia, Ariadne and 
Herwig, and the LUBOEI Bose-Einstein correlation model [21] implemented in Pythia, as 
discussed in Section 8.3. The effects of so-called 0{a) photon radiation have been studied 
using the KandY generator scheme [22], which uses YFSWW3 [23] and KoralW 1.51 [22] 
running concurrently, as discussed in detail in Section 8.4. 

All Monte Carlo samples have been passed through a complete simulation of the OPAL 
detector [24] and the same reconstruction and analysis algorithms as the real data. Small 
corrections were applied to the reconstructed jet and lepton four-vectors in Monte Carlo 
events better to model the energy scales and resolutions seen in data, as discussed in detail 
in Section 8.1. 

3 Event selection 

The selections of W^W" — qq^z^ and W"*'W~ — qqqq events are based on multivariate 
relative likelihood discriminants, and are discussed in detail in [25]. Events selected by the 
W'''W~ ivlu selection of [25] are rejected, and events selected as both qqiu and qqqq 
candidates are retained only for the qqiiy analysis. The sets of reference histograms used in 
the selections have been extended to maintain optimal performance for the highest energy 
LEP2 running. 

Scmilcptonic W^W~ qqiu decays comprise 44% of the total WW cross-section, and 
are selected using separate likelihood discriminants for the qqez^, qq/iv and qqri^ channels. 
These events are characterised by two well-separated hadronic jets, large missing momentum 
due to the escaping neutrino from the leptonic W decay, and in the case of qqez^ and qqnv 
decays, an isolated high-momentum charged lepton. In W+W~ qqTiy events, the r-lepton 
is identified as an isolated low multiplicity jet, typically containing one or three tracks. A 
small number of 'trackless-lepton' qqeu and qq/xz^ events are also selected, where the lepton 
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is identified based on calorimeter and muon chamber information only, withoTit an associated 
track. These events make up 2.5% of the qqez^ and 4.7% of the qq/xz^ samples. Hadronic 
W"'"W~ — > qqqq decays comprise 46 % of the total WW cross-section, and are characterised 
by four energetic hadronic jets and little or no missing energy. The dominant background 
results from e"^e^ — ^ Z/7 — ^ qq events giving a four-jet topology (qq — >■ qqqq or qqgg), and 
this is largely rejected using an event weight based on the 0{a1) QCD matrix element for 
this background process. 

The number of events selected in each of the channels and data-taking years is given 
in Table 1, together with the expectation from the Monte Carlo simulation with the WW 
production cross-section scaled to the prediction of KandY (which is more accurate than 
that of KoralW). The average selection efficiency and purity of each channel in the desired 
WW signal topology are also given, estimated from Monte Carlo events and averaged over 
all centre-of-mass energies. The dominant backgrounds are events misclassified between the 
qqeiy / qqiiu and qqriy channels in the qqii' selection, and e'''e~ ^ Z/7 ^ qq events giving a 
four-jet topology in the qqqq channel. Combining all three qqiv sub-channels, and including 
events mis-classified between them, 87 % of qq£u events are selected for the mass and width 
analyses. However, not all selected events are actually used by each analysis — some poorly 
reconstructed events are removed by analysis-specific cuts as discussed below. 

4 W boson reconstruction and kinematic fitting 

All three analysis methods use similar event reconstruction and kinematic fit techniques to 
determine the W mass on an event by event basis. In W^W~ qqli' events, the procedure 
begins by removing the tracks, electromagnetic and hadronic calorimeter clusters correspond- 
ing to the lepton identified by the event selection. A matching algorithm is then applied to 
tracks and calorimeter clusters, and the cluster energies are adjusted both to compensate for 
the expected energy sharing between the electromagnetic and hadronic calorimeters, and to 
account for the expected energy deposits from any associated tracks. This procedure has been 
optimised to obtain the best possible jet energy resolution on Z ^ qq events at ^/s « 91 GeV, 
where use of the hadronic as well as the electromagnetic calorimeter information improves 
the energy resolution by about 10%. The reconstructed objects (referred to hereafter as par- 
ticles) are then grouped into two jets using the Durham jet-finding algorithm [26]. Estimates 
of the jet energies, directions and masses are derived from the four-momentum sum of all the 
tracks and corrected calorimeter clusters assigned to the jet, assigning tracks the pion mass 
and clusters zero mass. Corresponding error matrices are also assigned to the reconstructed 
jet energies and directions, based on studies of jet resolution in Monte Carlo. 

In qqez^ events, the electron energy is reconstructed from the energy of the associated 
electromagnetic calorimeter cluster, and the direction is taken from that of the associated 
track (except in trackless qqeu events, where both the energy and direction are taken from the 
calorimeter cluster). In qq/v,zv events, the track is used for both the muon energy and direction 
estimates. In both cases, calorimeter clusters which are not associated to the lepton, but which 
are close to the lepton track and consistent with originating from final-state radiation, are 
added into the lepton energy estimate. In qqri^ events, the r energy cannot be reconstructed 
due to the undetected neutrino(s) produced in its leptonic or hadronic decay. This means 
that only the hadronic W ^ qq decay carries usable information about the W mass, and the 
r energy and direction are not reconstructed; this is also the case for trackless qqpLu events. 
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However, a complication can arise in the case of hadronic r decays if the r decay products 
are incorrectly identified and some of them mistakenly included in the reconstruction of the 
qq system. The mass information in such events can sometimes be recovered by using an 
alternative r reconstruction, forcing the whole event to a three-jet topology and assuming 
the r to be the jet with lowest invariant mass. A multivariate procedure based on angular 
and momentum variables is therefore used in hadronic r decays to decide between the two 
alternatively reconstructed topologies. 

In W+W~ qqqq events, the initial reconstruction used in the event selection is made 
by grouping all tracks and clusters into four jets using the Durham algorithm, with double- 
counting corrected as discussed above. However, a hard gluon is radiated from one of the 
quarks in a significant fraction of qqqq events, and the mass resolution for such events can 
be improved by reconstructing them with five jets [7]. The convolution and reweighting fits 
treat all qqqq events in this way, whilst the Breit-Wigner fit reconstructs the event as four 
or five jets depending on the value of ^45, the value of the Durham jet resolution parameter 
at which the five- to four-jet transition occurs. In all cases, the jets can be assigned to the 
two W bosons in several possible ways, leading to combinatorial background where the wrong 
assignment has been chosen — this is dealt with in different ways by the different analysis 
methods as discussed in detail below. 

The invariant masses of the two W bosons in the event could be determined directly 
from the momenta of the reconstructed jets and leptons, but the resolution would be severely 
limited by the relatively poor jet energy resolution of aEj^^/Ejet ~ 12% for well-contained 
light-flavour jets. For events without significant initial-state radiation, the W mass resolution 
can be significantly improved by using a kinematic fit imposing the four constraints that the 
total energy must be equal to the LEP centre-of-mass energy and that the three components 
of the total momentum must be zero (referred to as the 4C fit). Since the uncertainty on the 
two reconstructed W boson masses is typically still larger than the intrinsic W boson width 
of around 2 GeV, the resolution can be further improved by constraining the two masses to a 
common value (the 5C fit) . The 4C and 5C fits are used in various ways by the three analysis 
methods. In the qqeu and qqpiv channels three of the constraints are effectively absorbed by 
the unmeasured neutrino, and in the qqrz^ channel an effective one-constraint fit is performed 
to the hadronic part of the event only. In all kinematic fits, the velocity of the jet (3 = pjetZ-Ejet 
is kept fixed as the jet energy Ej^t is varied, which results in the jet momentum pjet and mass, 
mjet = \J -Ejet ~ pfet ' varying. This procedure is found to give results which are about 
1 % more precise than the fixed mjet approach used previously [7] . 

The dominant systematic error on the measurement of the W mass and width in the 
qqqq channel comes from possible final-state interactions between the decay products of the 
two W bosons. According to phenomenological models, these interactions mainly affect low 
momentum particles produced far from the cores of the jets. The uncertainties due to final- 
state interactions can therefore be reduced by deweighting such particles when calculating the 
jet four-momenta, for example by removing all particles with momentum p below a certain 
cut, weighting particles according to their momentum or only using particles whose directions 
lie close to the jet axis. 

Such an approach is used for the qqqq channel W mass measurement in this paper. The 
jet energy and mass are calculated using the original Durham jet definition, but the jet 
direction is taken instead from the sum of the momenta of all particles assigned to the jet 
which have p > 2.5 GeV. This cut strongly reduces the systematic uncertainties due to final- 
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state interactions, at the expense of some loss of statistical precision due to the reduction 
in jet angular resolution. This value of the cut was found to be optimal given the expected 
statistical error of the OPAL analysis. In around 4 % of jets, no particles have momenta above 
2.5 GeV, in which case the original jet direction is used. For comparison, the qqqq analysis 
results are also given using the unmodified Durham jet direction reconstruction (referred to 
as Jo), though this value is not used in the final result. In the qq£z^ channel and for the W 
width analysis, the unmodified Jq Durham jet reconstruction is always used. 

The sensitivity of the qqqq W mass analysis to final-state interactions can also be in- 
creased, by using a jet direction reconstruction giving higher weight to soft particles. In the 
convolution analysis, this is done by using a second modified reconstruction method, where 
the jet direction is calculated from the vector sum of the momenta of all particles assigned 
to the jet, each one weighted by p'^, with k = —0.5. The difference between the W mass 
calculated using this algorithm (referred to as k-o.s) and the algorithm with p > 2.5 GeV 
(referred to as P2.5) is sensitive to the presence of final-state interactions, and is used to set 
a limit on their possible strength within specific models. Using the same method, but with 
positive values of k, reduces the sensitivity of the analysis to final-state interactions, as does 
using a cone-based direction reconstruction where only particles within an angle R of the 
original jet axis are used to calculate an updated jet direction. Results from these algorithms 
are also given in Section 5.3 for comparison purposes. 



5 The convolution fit 

The convolution fit is based on the event-by-cvcnt convolution of a resolution function i?j(m/^, 7712) 
for event i with a physics function P{m'i,m'2\m-w, Fw, V^)- The latter represents the expected 
distribution of true event-by-event W masses m'^ and m2 given the true W mass and width, 
and the centre-of-mass energy. The resolution function gives the relative probability that a 
given observed event configuration could have arisen from an event with true masses m'^ and 
and is calculated in different ways for the qqiv and qqqq channels. The physics function 
is the same for both channels, and is given by 

P(m'i, 771-2 I "^w, Tw, = (1) 



Oo 



■B{m'i I 777W, Tw) • B{m'2 \ 777w, Tw) • S{m[,m'2 \ Vs') (8) I{y/s, \fs' 



where ao normalises the integral of P over the (777'^, 7772) plane to unity and the symbol '(8)' 
denotes convolution. The unnormalised relativistic Breit-Wigner distribution B is given by 

777^ 

5(777 I Tnw, rw) = (^2 _ ^^2)2 + (^2r^/^^)2 ' (2) 

and the phase space term S, describing the suppression close to the kinematic limit 777^+7772 = 
\fs', where \/s' is the effective centre-of-mass energy after initial-state radiation, is given by 

5(777;, 777'2 I = ^J{S' - (777; + 777'2)2) • (s' - (777^ - 777'2)2) . (3) 



The radiator function I{^/s, v s') describes the effect of initial-state radiation causing an event 
of centre-of-mass energy ^/s to have its effective centre-of-mass energy reduced to ^/s' and is 
given by 



IiV-s,V?)=Px^^-^^-^^^^ (4) 
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where (j{-\fs,my^) is the W-pair production cross section for a given y/s and mw, x is 
the normahsed initial state radiation photon energy x = E^/^/s, = (s — s')/2y^, and 
P = {2a/'K)\og{{y/s/m^y' — 1) where a is the electromagnetic coupling constant and the 
electron mass [10]. 

The signal likelihood for event z, >Cf^(mw,rw)) is calculated from the convolution of the 
resolution and physics functions: 

Cf^{mw, Tw) = Ri{m[,m2) P{m[,m2 \ mw, Tw, -v/sj ■ (5) 

Additional terms Cf^ and cf^"' are included to account for the presence of background from 
ZZ and Z/7 production. These likelihoods are parameterised using Monte Carlo events and 
weighted by event-by-cvent probabilities V^^^, Vf^ and "Pf^^ that the event comes from each 
of these sources, derived from the event selection likelihoods. The total likelihood for event i 
is then given by 

£r (mw, Tw) = T^r-Crimw, Fw) + V^Cf + vf'^C^J^ , (6) 

and the likelihood for the whole sample is given simply by the product of the individual event 
likelihoods. The convolution integrals in Equations 1 and 5 are performed numerically, and 
evaluated using a grid of 8100 points in the part of the (777,1,777,2) plane satisfying 100 GeV < 
m'l + 77^2 < y/s and |777^ — 777,2! < GeV. 

Separate fits axe performed to extract the W mass and width. For the mass, 777w is varied 
to maximise the overall likelihood, with Fw determined from 777w by the Standard Model 
relation [10] 

Fw = 3Gf 777^(1 + 2a,/37r)/(2V2^) , (7) 

where Gf and ag are the Fermi and strong coupling constants. The fitted mass is obtained 
from the maximum of the likelihood curve, and then corrected for the biases discussed below. 
For the W width, th-w is kept fixed at 80.33 GeV and only Fw is varied. In the qqqq channel, 
the fitted mass does not depend on the assumed width and vice versa, but in the qq^T^ channel 
the width has a small residual dependence on the assumed W mass. This is corrected at the 
end of the fit procedure according to the value derived from the mass fit, a simultaneous 
two-dimensional fit of m-w and Fw not being possible for computational reasons. 

5.1 The qq£i/ convolution fit 

In W"'"W^ qqiv events, the missing neutrino leads to kinematic fit solutions with likeli- 
hoods which are not Gaussian, especially if the constraint that the two W masses are equal is 
not applied. The convolution fit provides a natural framework to exploit all available infor- 
mation in the non-Gaussian resolution function i?(777^ , 7772)- For each event, this function is 
mapped out in the (777/^ , 777^2 ) plane by performing many six-constraint kinematic fits, where in 
addition to energy-momentum conservation, the two W masses are fixed to the input values 
m'l and 7773 rather than being left free to be determined in the fit. Each fit therefore gives 
only a value, which varies as a function of m'l and 7772 and expresses the consistency of the 
event with the input W mass hypothesis. The minimum (x^j^) of this contour corresponds 
to the fitted values of the two W masses, 7711 and 7712, which would have been returned by 
a standard 4C fit. The resolution function R at each point is derived from the contour 
via the relation R = e^UXmin ~ X^)/^) and normalised so that its integral is unity over the 
{m'l, 7712) plane. 
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The kinematic fits in the qqei^ and qq/Ui' channels are performed using semi-analytic 
approximations with two simplifying assumptions, namely that the lepton direction is fixed, 
and that the fitted jet directions are constrained to lie in the plane defined by their measured 
values. These allow the fit to be reduced to a one-dimensional numerical minimisation. If 
the event is very badly measured (or is in fact a background event) , the apparent minimum 
of the contour may lie at the edge of the mass grid — this happens in about 5 % of qqezv 
and qq/iv candidates. An attempt is made to recover some useful W mass information from 
such events by discarding the lepton and refitting them as qqru events (where the lepton 
information is never used); this is also done for trackless-lepton qqnf candidates which have 
no useful estimate of the lepton energy and only a poor estimate of its direction from the 
muon chambers. 

The kinematic fit for qqru events involves only the hadronic system, the only variables 
of interest being the angle between the two jets from the W — >^ qq decay and the sharing 
of the available beam energy between them. The resolution function R is mapped out using 
the same technique as for qqev and qqjiv events. Events with a solution within 2.5 GeV of 
any edge of the mass grid are not considered further; this happens to about 25 % of signal 
W+W~ qqriy events and 58 % of background qqri^ candidates. 

The non-WW sources of background in the qq^z^ channel are very small in all but the 
qqr:^ case. Their contributions are accounted for by background terms in the likelihood (see 
Equation 6) which are parameterised as functions of mi and m2, the two W masses at the 

minimum of the kinematic fit solutions, and y^. Separate parameterisations are used for 
qqez^, qq/ii^, qqriy and trackless-lepton qqcz^ and qq/J^v candidates, derived from large samples 
of background Monte Carlo simulated events. 

In Monte Carlo qq£u events, the W mass and width estimates derived from the convolution 
fit differ from the simulated values by up to 350 MeV, due to effects not fully accounted for 
in the likelihood function, for example biases in the input jet and lepton four-vectors, and 
imperfections in the treatment of initial-state radiation and backgrounds. These biases are 
studied by applying the convolution fit to large Monte Carlo samples of simulated signal and 
background events with various true values of mw from 79.33 GeV to 81.33 GeV, and Fw 
from 1.6 GeV to 2.6 GeV. For the W mass fit, the biases are found to depend on ^/s but 
not on the true values mw and Fw, and are parameterised from Monte Carlo as smooth 
functions of ^/s. In the width fit, the bias on the reconstructed width is found to depend 
slightly on the true width, as well as on the true mass as discussed above. These biases are 
again parameterised using Monte Carlo. The errors returned by the fits are also checked, by 
studying pull distributions obtained from fits to many Monte Carlo subsamples constructed 
so as to have the same integrated luminosity as the data in each year. These studies show 
that the fits underestimate the statistical error by about 5 %, reflecting imperfections in the 
input jet and lepton error matrices. Corresponding corrections are therefore applied to the 
statistical errors determined by the fits. 

After these corrections, the fits give unbiased results on Monte Carlo samples, but several 
further small corrections, amounting to a total of about 5 MeV for the W mass and 20 MeV 
for the W width, are applied to the data results. These account for effects not present in the 
Monte Carlo samples used to calculate the bias corrections, namely additional non-simulated 
detector occupancy, deficiencies in the description of kaon and baryon production in the JT 
hadronisation model, and 0{a) photon radiation modelled by KandY but not KoralW. 
These corrections, which are also applied to the results of the reweighting and Breit-Wigner 
fits, are discussed individually in more detail in Section 8. 
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A second convolution-based fit method (referred to as tfie 'CV5' fit) is used to make an 
additional cross-cfieck of the W width fit result in the qqiu channel. This method is similar to 
the fit described above, except that the two input W masses m'l and m2 are set to be equal, 
making it equivalent to a 5C rather than a 4C fit, and tracing out the probability contour 
only along the diagonal m'l = m'2 in the (m'^^, 77X2) plane. This reduces the number of kinematic 
fits needed per event, allowing the values to be determined using numerical minimisation 
rather than the fast analytic approximations used above. In this fit, a single Breit-Wigner 
distribution is used in the physics function analogous to Equation 1, and both the mass 
mw and width Tw are determined simultaneously. Similar bias correction procedures and 
parameterisations are used as for the standard convolution fit. 

5.2 The qqqq convolution fit 

The qqqq channel differs from the qq£i^ channel in several important respects: no prompt 
neutrinos are produced, leading to better constrained kinematics, but the assignment of jets 
to the two decaying W bosons is ambiguous, leading to combinatorial background where the 
wrong assignment is made. Non-WW background (particularly from e+e" 7j/"f events 
producing four jets) is also much more important than in qqiv events, contributing 16 % of 
the selected sample. 

In a significant fraction of W"^ W~ qqqq events, a hard gluon is radiated from one of the 

quarks, and these events are better reconstructed as five-jet rather than four-jet events. Since 
the division between four and five jets is rather arbitrary, the convolution fit reconstructs all 
qqqq events with five jets. In qqqq events with no hard gluon radiation one of the quark 
jets is split in two by this procedure, but the two jet fragments have a high probability to 
be correctly assigned to the same W boson. A more serious problem is the combinatorial 
background — with five jets there are ten possible assignments of the jets to two W bosons, 
compared with only three in a four-jet topology. This is dealt with in two ways. Firstly, only 
4C fits are used, where the two W boson masses are not constrained to be equal; many of 
the incorrect jet assignments give kinematic fit solutions with two very different masses, in 
contrast to the correct solution with two similar masses. Secondly, energy ordering of the 
jets is used together with an artificial neural network algorithm based on the 4C fit mass 
differences to weight each remaining jet assignment combination in the likelihood fit. 

In more detail, an initial 4C kinematic fit imposing four-momentum conservation is applied 
to the five jets, which are then ordered according to their fitted energies. The event is also 
reconstructed in a four-jet topology using the Durham scheme, resulting in two of the five 
jets being combined, the other three remaining unchanged. The three unchanged jets are 
labelled 1-3 such that Ei > E2 > E^, whilst the remaining two jets in the five-jet topology 
arc labelled 4 and 5, with E/i > E^, where Ei refers to the fitted energy of jet i from the 
initial 4C kinematic fit. These jets can be assigned to two W bosons in ten different ways, 
with the combinations numbered (124,35), (125,34), (12,345) and so on. The combination 
(123,45) is not considered further, as it has one W boson formed from just the split jets, 
which is very unlikely. For each of the nine remaining combinations c, the jet four-vectors 
resulting from the 4C fit are combined to calculate the reconstructed masses and 777,2 ' and 
the mass difference drric = 777f — 7772- The nine mass differences are input to an artificial neural 
network [27] with seven outputs, corresponding to each of the remaining combinations apart 
from (124,35) and (134,25), which are also discarded at this stage, having little probability of 
being correct due to the large imbalance in the energies assigned to the two W bosons. The 
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network is trained using a large sample of signal WW Monte Carlo events to give values close 
to one at the output corresponding to the correct combination, and zero for all other outputs. 
In each qqqq event, the seven outputs are normalised to sum to unity, and all combinations 
with output Qc > 0.12 are retained for the final likelihood fit. This cut value is found to 
minimise the statistical error on the W mass in Monte Carlo events. 

The distribution of Qc for all jet combinations with Qc > 0.05 in data is shown in Fig- 
ure 1(a), together with the expectation from Monte Carlo, broken down into correct and 
wrong combinations in WW events, and ZZ and Z/7 background. The fraction of WW 
Monte Carlo jet combinations which are correct^ is shown as a function of Qc in Figure 1 (b) . 
The probability for each of the jet combinations to be correct before the neural network se- 
lection is shown in Figure 1(c), showing the power of the initial energy ordering in already 
distinguishing the correct combination. The number of combinations A'comb with Qc > 0.12 
is shown in Figure 1(d) — typically 3 or 4 combinations are retained for the final fit. Some 
discrepancies between data and Monto Carlo are visible; these are addressed in the systematic 
uncertainty studies as discussed in Section 8.7. 

The resolution function for each retained combination is generated from the fitted W 
masses mi and m2 returned by the 4C kinematic fit, together with their associated errors 
and correlation coefficient. This simple approach, rather than mapping out the full resolution 
function using many six-constraint fits with fixed input W masses, is adequate due to the 
better-constrained kinematics compared with the qqiu channel. However, to model the tails 
better, a two-dimensional double Gaussian resolution function is used, with separate core and 
tail components. The core has a width given by the evcnt-by-cvent kinematic fit errors and 
a weight of 58 %, whilst the tail component contributes the remaining 42 % of the resolution 
function and has a width 2.2 times larger than the core. These associated parameters were 
derived from studies of the fit resolution in Monte Carlo simulation. 

The signal likelihood function is more complicated than that for qqii' events as it must 
account for the several jet assignment combinations in each event. This is achieved by treating 
the W boson masses mf and for each combination c as independent observables. For each 
combination, one pair of masses is described by the convolution of signal resolution and 
physics functions, whilst the others (considered to come from combinatorial background) 
are each described by a parameterised function C(mi, m2, y^), where (j+ is the error 
on the sum of the two masses mi -|- m2- This function is obtained from the distributions 
of combinatorial background combinations in Monte Carlo events. The likelihoods for each 
of the iVcomb combinations are then summed, weighting each one by its associated neural 
network output Qc- Thus, the signal likelihood for one event is given by 

(mw,rw) = E Qcdc (8) 

c=l 

-^comb 

The overall event likelihood is again given by Equation 6. In the case of the qqqq channel, 
the likelihood CT'^ for ZZ events is also given by Equation 8, with mw replaced by the 

■'A jet combination is considered to be correct if all the jets are assigned to the correct W bosons, a jet being 
correctly assigned if more than half of its energy results from the decay products of the associated boson. In 
about 2 % of selected events no combination is considered correct, due to more than three jets being assigned 
to one boson according to this definition. 
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Figure 1: Details of jet assignment for the qqqq convolution analysis: (a) distribution of 
neural network outputs Qc for jet combinations in data and Monte Carlo, showing the correct 

and wrong combinations in WW events (separated by the dashed line), the contributions 
from ZZ and Z/7 background events, and the cut used to select combinations for fitting; 
(b) fraction of correct combinations in Monte Carlo WW events as a function of Qc, (c) 
probabilities that each jet combination is correct, based on energy-ordering before the neural 
network selection (note that combinations (124,35) and (134,25) arc never used in the fit); 
(d) number of accepted combinations A'comb per event in data and Monte Carlo. 
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Reweightin 
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Breit-Wigner 


Channel 


Fitted mw 


•^cxp 


Fitted mw 




Fitted mw 






(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


qqei^ 


80.511 ± 0.084 


0.088 


80.492 ± 0.088 


0.091 


80.500 ± 0.100 


0.104 


qqfiu 


80.432 ± 0.090 


0.089 


80.488 ± 0.093 


0.091 


80.523 ±0.105 


0.106 


qqriy 


80.354 ±0.126 


0.125 


80.289 ±0.130 


0.130 


80.269 ±0.135 


0.140 


qqiu 


80.451 ± 0.056 


0.056 


80.451 ± 0.057 


0.058 


80.457 ± 0.064 


0.065 


qqqq fe.s) 


80.353 ± 0.060 


0.059 


80.308 ± 0.064 


0.066 


80.286 ± 0.073 


0.075 


qqqq (Jq) 


80.394 ± 0.051 


0.051 


80.383 ± 0.056 


0.057 


80.424 ± 0.059 


0.060 


qqqq (k-o.s) 


80.508 ± 0.073 


0.073 











Table 2: W mass results (with statistical errors only) for each channel and fitting method. 
The expected statistical errors o"cxp from Monte Carlo subsample tests are also given. In the 
qqqq channel, results are given for the P2.5 jet direction reconstruction, the Jq reconstruction 
and the K-0.5 reconstruction giving increased weight to low-momentum particles (the latter 
for the convolution fit only). 



known mz, and the likelihood C'^^"' for Z/7 events is given by a similar expression but with 
no correct combination, only terms involving combinatorial background C{mi,m2,cr+, y/s). 
The parameterised functions C are determined from Monte Carlo separately for each type of 
event, and the event type probabilities V^^^, V^^ and T^^/f are also parameterised, as linear 
functions of the event selection likelihood. 

As for the qq^z^ fit, bias corrections are applied to the raw mass fit results, parameterised 
as a function of ^/s. These corrections are calculated separately for the fits using the Jq and 
modified jet direction reconstruction methods, and are largest (up to 400 MeV) for the K-0.5 
jet reconstruction. No significant dependence of the corrections on the true value of mw 
is observed. In the qqqq channel, the width fit bias is also found to be independent of the 
true value of Fw, and on the assumed value of mw- Monte Carlo subsample tests are also 
performed, and small corrections to the fit error estimates of typically 5-10 % are derived. 
Further small corrections of up to 9 MeV arc applied for effects not present in the default 
Monte Carlo samples, as discussed in Section 5.1. 

5.3 Convolution fit results 

The convolution fit is used to analyse the data for each year separately, and the results are 
then combined. The results and associated statistical uncertainties are given in Table 2, for 
the qqeu, qq/Ui/, qqrz/, combined qqiu and qqqq channels. In the qqqq channel, the results 
are given for jet direction reconstruction methods P2.5 and k_o.5 and for comparison also 
with the unmodified Durham jet algorithm (Jq) as used in the qq^z^ channel. The quoted 
results include all corrections made to the fit results as discussed above, but the averages do 
not include the effects of systematic uncertainties (the final results including all uncertainties 
are given in Section 9.2). Table 2 also gives the expected statistical errors for each channel, 
evaluated using fits to many Monte Carlo subsamples, each constructed to have the same 
integrated luminosity and centre-of-mass energy distribution as the data. In all cases, the 
data statistical errors are consistent with the expectations from Monte Carlo, after taking 
into account the expected level of statistical fluctuations. 

Distributions of the mean of the two W masses reconstructed in each event ((mi ±m2)/2) 
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Channel 


Convolution 

Fitted Fw Ccxp 
(GeV) (GeV) 


Reweighting 

Fitted Fw fcxp 
(GeV) (GeV) 


CVS convolution 

Fitted Fw CTcxp 
(GeV) (GeV) 


qqei' 

qqrz/ 


1.696 ±0.202 0.204 
2.181 ±0.233 0.216 
1.763 ± 0.289 0.289 


2.009 ± 0.200 0.202 

2.146 ±0.224 0.227 
2.089 ± 0.276 0.309 


1.975 ± 0.230 0.224 

2.138 ± 0.233 0.218 
2.204 ±0.188 0.241 


qqiu 

qqqq (Jq) 


1.926 ±0.135 0.134 
2.125 ±0.111 0.114 


2.088 ±0.131 0.127 
2.176 ±0.129 0.134 


2.103 ±0.120 0.131 



Table 3: W width results (with statistical errors only) for each channel and fitting method. 
The expected statistical errors Uexp from Monte Carlo subsample tests are also given. The 
Breit-Wigner fit does not measure the W width. 



are shown for the qqiv fit and sub-channels in Figure 2, and for the qqqq fit (with jet direction 
reconstruction P2.5) in Figure 3. In the latter figure, reconstructed mass combinations are 
shown for two ranges of Qc, showing the suppression of the combinatorial background achieved 
by the neural network algorithm. The results obtained in each year of data-taking are shown 
as the 'CV points in Figure 4; all the results arc consistent with the overall mean for each 
channel, and the values for the qqiv and qqqq {P2.5) averages are 4.3 and 1.0, each for four 
degrees of freedom. 

The corresponding results for the width are shown in Table 3 and Figure 5. The 1996 data 
at i/i ~ 172 GcV arc not used for the width analysis. Again, the statistical uncertainties are 
compatible with expectations, and the individual year results are consistent, the values 
being 5.4 and 2.1 for the qqiv and qqqq channels, each for three degrees of freedom. The 
statistical correlation between the qqlu mass and width results is estimated using Monte 
Carlo subsamples to be —0.19, whereas that for the qqqq channel results is found to be 
negligible. Note that the qqqq width analysis is performed using the unmodified Jq jet 
direction reconstruction as this gives the optimal balance between statistical and systematic 
errors from hadronisation and final-state interactions, and minimises the total error. The 
width result from the CVS convolution fit in the qq^z^ channel is also shown; this fit also 
measures the W mass and gives a result of 80.424 ± 0.077 GeV, consistent with that derived 
from the standard convolution fit. The statistical correlation coefficient between the CVS 
width and mass fit results is 0.28. 

In order to study the evolution of the fitted W mass with changing jet direction recon- 
struction, the complete convolution fit has been repeated fifteen times, using momentum cuts 
at 1.0, 1.75, 2.0, 2.5, 3.0 and 4.0 GeV (2.5 GeV being used for the qqqq analysis result in this 
paper), momentum weights p'^ with k values of 0.5, 0.75, 1.0, —0.5 and —0.75, and cones of 
half-angle R =0.3, 0.4, 0.5 and 0.6 rad. For each jet direction reconstruction method, the 
mass difference with respect to the Jq direction reconstruction using all particles associated 
to the jet Am{Jx, Jo) = m{X) — m(Jo) is calculated in both the qqiv and qqqq channels. 
The statistical error on Am{Jx, Jo) is also calculated, using Monte Carlo subsamples to take 
into account the correlation due to the different reconstruction methods being applied to the 
same events. The results are shown in Figure 6. The qqiu Am{ Jx , Jo) values are generally 
slightly positive, but the changes in fitted W mass are consistent with the expected level of 
statistical fluctuations, demonstrating that the qqiu W mass results are stable with respect 
to changing the jet direction reconstruction over all methods and a wide range of parameter 
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Figure 2: Reconstructed mean mass distributions for W+W" — > qqiv, qqev, qq/iv and qqriy 
candidates fitted using the convolution analysis. The points with error bars show the data, 
and the histograms show the Monte Carlo expectation (with mw = 80.415 GcV), broken 
down into contributions from signal WW events with the correct Icpton type, WW events 
with mis-identified leptons ('WW mis-ID'), and background from ZZ and Z/7 events. 
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Figure 3: Reconstructed mean mass distributions for W^W^ qqqq combinations fitted 
using the convolution analysis with the P2.5 jet direction reconstruction, for two different 
ranges of correct combination probability Qc- The points with error bars show the data, 
and the histograms show the Monte Carlo expectation (with mw = 80.415 GeV), broken 
down into contributions from correct jet combinations in signal WW events, combinatorial 
background in WW events, and background from ZZ and Z/7 events. 



values. The qqqq Am(Jx, Jo) results also tend to be close to zero, with the exception of the 
result from jet direction reconstruction k_o.5 (with enhanced sensitivity to low momentum 
particles), which is significantly higher (Am(K_o.5) <^o) = 114ib47MeV) than the results with 
all other jet direction definitions. The result from method a-lso shows a high value 

of Am(«;_o.75, Jo); although with low significance. Note that hadronisation uncertainties are 
also significant in these comparisons, and increase to around 20 MeV for the alternative jet 
direction reconstruction methods, as discussed in Section 8.2. 

Sensitivity to final-state interactions in the qqqq channel can be maximised by studying 
the variable 

Am{Jx, K-0.5) = m{X) - m(/c_o.5) , (9) 

the difference in W mass between jet direction reconstruction methods with reduced and in- 
creased sensitivity to these effects. The largest deviation from zero is seen for Am(p2.5, 1^-0.5) = 
— 152 lb 68 MeV, where the error is purely statistical, but takes into account correlations be- 
tween the different reconstruction methods. The use of the mass differences to place limits 
on the effect of final-state interactions, specifically colour reconnection, is discussed in Sec- 
tion 9.1. 
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Figure 4: W mass results (with statistical errors only) for each channel in the convolution 
(CV), reweighting (RW) and Breit-Wigner (BW) fitting methods, as a function of data- 
taking year, and for all years combined. Results are shown for the qqei^, qq/iz^, qqriy and 
combined qq^z^ channels, and for the qqqq channel with the modified (P2.5) and Jq jet direction 
reconstruction methods. The results for the Breit-Wigner qqlv and qqqq fits to 1996 data are 
taken from [5] . The dotted line indicates the central measured value of mw from the present 
paper. 
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Figure 5: W width results (with statistical errors only) for each channel in the convolution 
(CV), reweighting (RW) and 5C convolution (CV5) fitting methods, as a function of data- 
taking year, and for all years combined. Results are shown for the qqei^, qq/iz^, qqriy and 
combined qqiv channels, and for the qqqq channel with the Jq jet algorithm. Only the 
reweighting fit measures the W width using the 1996 data, where results are shown for the 
qqiu and qqqq channels only. The dotted line indicates the central measured value of Fw 
from the present paper. 
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Figure 6: Differences in mw measured using various jet direction reconstruction methods 
X and the Jq reconstruction method, for both qq£zv and qqqq data in the convolution fit. 
The points are highly correlated; the errors are purely statistical and take into account the 
correlation between the result from each alternative direction reconstruction method and that 
from Jq. 



6 The reweighting fit 

The reweighting fit extracts the W mass and width by comparing reconstructed data dis- 
tributions with Monte Carlo 'template' distributions with varying mw and Fw- Templates 
of arbitrary mw and Tw are obtained by reweighting Monte Carlo simulated data samples 
containing all signal and background final states, and a maximum likelihood fit is used to 
find the values of mw and Fw that best describe the data. The reweighting fit is a more 
sophisticated development of that used in [7], the main changes being the use of simultaneous 
reweighting in three (qqez^, c^qpLV and qqqq) or two (qqrz^) reconstructed variables, and an 
improved procedure for handling the combinatorial background in the qqqq channel. 
In more detail, the likelihood for each event i is given by 

CT\ai\m^, Fw) = ^^^^^^^(^.i^^^ r^) + pzz^zz^^.^ ^ pZ/^^z/^^^.) ^ (^q) 

where a is the set of reconstructed variables used in the likelihood, £'^^s(a|mw) Tw), C^'^{a) 
and C'^'^{a) are the likelihood distributions of the variables a in WW, ZZ and Z/7 events, 
and r^'^, r'^'^ and V'^'^ are the (fixed) fractions of WW, ZZ and Z/7 events in the sample, 
estimated from Monte Carlo simulation at the corresponding centre-of-mass energy. The 
signal probability distribution >C^'S(a|mw, Tw) is obtained by reweighting four-fermion Monte 
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Carlo events with a true W mass of = 80.33 GeV and a width of = 2.09 GeV by the 
ratio of two Breit-Wigner functions. The weight /, for a Monte Carlo event i with true 
event-by-event W boson masses m'l and m'2 is given by 



where the Breit-Wigner function B{m\m^J^^T^J^l) is given by Equation 2. The effect of back- 
ground is accounted for via the background terms in Equation 10, whose probability distribu- 
tions are calculated as a function of a. using large samples of unweighted background Monte 
Carlo events. 

The probabihty distributions C^^^, and C^^'^ are calculated in bins of the reweighting 
fit variables a, the bin size varying with a in order to achieve an approximately constant 
number of events per bin and minimise fluctuations from limited Monte Carlo statistics. The 
likelihood for the whole sample is therefore obtained from the number of events Nj in each 
bin j where the variables take the values df 



Two types of reweighting fit are performed. In the first, the likelihood >C*°* is maximised as 
mw is varied and Fw is determined from mw by the Standard Model relation (Equation 7). 
In the second, a two parameter fit is performed, allowing both mw and Fw to vary simul- 
taneously. The results for mw are very similar in both cases, but for consistency with the 
convolution and Breit-Wigner analyses, the result of the first fit is used for the reweighting 
fit W-mass result in this paper. 

6.1 The qq^i/ reweighting fit 

The qq^z^ reweighting fit uses the same basic event selection as the convolution fit. In the qqez^ 
and (\qpLu channels, both 4C and 5C kinematic fits are performed, and all events for which 
both kinematic fits converge are retained. The reweighting fit is performed simultaneously in 
three reconstructed variables which make up the variable set a: 

• The reconstructed W mass from the 5C fit, msc, in 16 bins from 65 to 105 GeV. 

• The error on the reconstructed 5C fit mass, (7^^,^, in five bins from 0.5 to 6.5 GeV. 

• The two-jet invariant mass from the 4C fit, in four bins from 40 to 140 GeV. 

The bin sizes vary, and are chosen such that each of the 320 bins in each channel is populated 
by about 400 Monte Carlo events. Events arc discarded if any of the variables fall outside the 
bin ranges. The use of the error on the 5C fit mass and the jet-jet invariant mass significantly 
improves the statistical precision of the fit as compared to the one-dimensional reweighting 
using m5c alone [7]. 

In the qqTiy channel, all information comes from the hadronic system and is extracted 
using an analytic implementation of the 5C fit. Two variables are used, namely the 5C fit 
mass (20 bins from 65 to 105 GeV) and its error (five bins from to 6.5 GeV). A variable bin 
size is again used, with around 1000 Monte Carlo events per bin. Events with a kinematic fit 
probability of less than 10~^ are removed. 



i?(m']^ I mw, Fw) -^(mglmwjFw) 
B{m[\m'^,T%) i?(mymO^, F^^) 



(11) 




J 



(12) 
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The reweighting fit technique should imphcitly correct for all effects which bias the recon- 
structed W mass, providing they are included in the Monte Carlo simulation used to generate 
the template distributions. Therefore, the effects of initial-state radiation, event selection and 
reconstruction biases are all included. This is checked using large Monte Carlo samples over 
the full range of centre-of-mass energies and true W masses from 79.33-81.33 GeV. The errors 
returned by the fits are similarly checked by studying the pull distributions in Monte Carlo 
subsamples, and found to be unbiased. 

6.2 The qqqq reweighting fit 

The qqqq reweighting fit uses the same basic event selection as the convolution fit. However, 
the method used for the assignment of jets to the two W bosons is rather different, with only 
one combination per event entering the final fit. The tracks and clusters of the event are first 
grouped into five jets using the Durham jet algorithm, and a 4C kinematic fit is performed. 
The value of the variable yij = EiEj{l — cos 6ij) calculated for each pair of jets, where Ei and 
Ej are the fitted energies of jets i and j and 6ij is the angle between them. The five fitted jets 
are assigned to the two W bosons requiring that the pair of jets with the lowest y^j is always 
kept together, both jets being assigned to the same W boson. Each of the other three jets is 
then assigned in turn to the same W boson as the paired jets. This results in three distinct 
jet assignment combinations, which are each fitted with a 5C kinematic fit. 

In the qqqq analysis with the p2.5 jet direction reconstruction method, the best of the 
three jet combinations is determined using the jet-pairing likelihood technique described in 
[7], with the two jets corresponding to the minimum y^j merged into a single jet. For each 
of the possible jet pairing assignments, three input variables are calculated and fed into a 
likelihood discriminant, and the combination with the largest output value is retained for the 
fit. The likelihood reference distributions are determined using large Monte Carlo samples, 
separately at each centre-of-mass energy. The input variables consist of the value of the CC03 
matrix element for W-pair production [25] , determined from the measured four- vectors of the 
reconstructed jets; the difference in reconstructed masses of the two W bosons, determined 
using the initial 4C kinematic fit; and the sum of the di-jet opening angles. The CC03 matrix 
element is averaged over three assumed W mass values from 80.1 to 80.6 GeV. This algorithm 
selects a jet assignment combination in every selected qqqq event, and is correct 72 % of the 
timc.^ For the qqqq analysis with the Jq jet algorithm, the jet angular resolution is such that 
the CC03 matrix element provides good discriminating power by itself, and the jet-pairing 
likelihood is not used. This algorithm selects the correct jet assignment in 74 % of cases. 

Having selected one jet assignment combination, the corresponding 4C and 5C kinematic 
fits are used to provide the reconstructed variables entering the reweighting fit likelihood. 
These variables arc: 

• The reconstructed W mass from the 5C fit, m5c, in 24 bins from 65 to 105 GeV. 

• The error on the 5C fit mass, CTj^gc, in 5 bins from zero to 5 GeV. 

• The difference of the two 4C fit masses, 5m, in 5 bins from —50 to 55 GeV. The mass 
difference is signed such that the W boson containing the jet with the highest energy 
before the kinematic fit contributes with a positive sign. 

"'in this case, each jet is Eissociated to the original quark closest to it in angle, and the jet assignment is 
considered correct if all quarks associated to jets assigned to one W boson do in fact originate from the decay 
of one boson. 
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As for the qqiiy fit, tlie bin sizes are chosen so that each bin is populated by around 400 Monte 
Carlo events. Events in which either of the fits fail, or in which any of the reconstructed 
variables fall outside the given range, are discarded. 

The fit method is checked using large Monte Carlo samples as for the qqtv reweighting 
fit. The errors returned by the fit are also checked by studying pull distributions in Monte 
Carlo subsample tests, and found to be unbiased. 

6.3 Reweighting fit results 

The reweighting fit is used to analyse the data from each year and channel separately, and the 
results for the different years are combined to give the values shown in Tables 2 and 3. The 
results from each year are also shown separately as the 'RW points in Figures 4 and 5. As 
discussed above the mass values are determined using a one parameter fit to m-yv only, and 
the width values are determined using a two parameter fit to mw and Tw- In the latter fits, 
the correlation coefficients between mw and Fw are 0.08 in the qqiu channel and 0.07 in the 
qqqq channel, and the mass results agree with those from the one parameter fits to within 
1 MeV. No separate Fw results are shown for the individual qqciy, qqfiv and qqri^ channels in 
1996 due to the small numbers of selected events, but the 1996 data are included in the overall 
averages. The expected statistical errors are also given in Tables 2 and 3, evaluated using 
Monte Carlo subsample tests. The statistical errors on the data results are again consistent 
with those expected from Monte Carlo, taking into account the expected level of statistical 
fluctuations. 

The reconstructed mass distributions from the 5C fits can be seen in Figures 7(a) and 
(b), for both the qqlu and qqqq channels (for the latter only the selected jet assignment 
combinations are shown). The rewcighted Monte Carlo template distributions corresponding 
to the fitted values of mw in each channel are also shown, including both signal (WW events 
with the correct jet assignment) and background contributions. The width of the qqiu mass 
peak is smaller than that from the convolution fit shown in Figure 2 because the latter displays 
the average of the two fitted W masses in each event. This average does not take into account 
the better resolution of the qq system mass compared to that of the iv system, information 
which is however included in the convolution fit itself. 

7 The Breit-Wigner fit 

The Breit-Wigner fit is based on a simple likelihood fit to the distribution of W boson masses 
reconstructed using a 5C kinematic fit in each event, and is very similar to that described 
in [5]. The main motivation for this analysis is to extract the W mass using a simple and 
transparent method, to act as a cross-check for the convolution and reweighting fits. The 
Breit-Wigner fit docs not measure the W width. 

The event selection and reconstruction are very similar to those of the convolution and 
reweighting fits. In the qqiv channel, only events with a 5C kinematic fit probability exceeding 
10~^ are used in the analysis. Events in each of the lepton sub-channels (qqei^, qqpLf and 
qqru) are treated separately, and the qqrz^ channel is further divided into events where the r 
decays leptonically or hadronically. In the qqqq channel, events are reconstructed as five jets 
if the Durham jet resolution parameter 7/45 > 0.0037 (about 23% of the events), and as four 
jets otherwise. In four-jet events, 5C kinematic fits are performed on all three possible jet 
pairings. The fit with the highest probability Pi is used if Pi > 0.003 for the P2.5 jet direction 
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Figure 7: Distributions of 5C fit masses for the qqiv and qqqq channels in the reweighting 
(a,b) and Breit-Wigner (c,d) fits. For the reweighting fits, the histograms show the reweighted 
template distributions corresponding to the fitted W mass values, and for the Breit-Wigner 
fits, the fitted functions are indicated by the solid lines drawn over the fit regions (70-88 GeV). 
In both cases, the estimated background contributions are indicated by the shaded regions. 
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reconstruction method, and Pi > 0.01 for the Jq method. The fit with the second-highest 
probabihty P2 is also used (with equal weight) if it passes both the previous probability cut 
and P2 > g-Pi; this occurs in approximately 20% of events. In five-jet events, at most one 
of the possible ten jet assignment combinations is used, selected according to the output of 
the jet assignment likelihood algorithm used in [7]. The likelihood inputs are the difference 
between the two W masses in a 4C fit, the largest inter-jet opening angle between jets in the 
three-jet system, and the cosine of the polar angle of the three-jet system. The jet combination 
giving the largest likelihood value is used provided the value is greater than a minimum cut 
requirement, which happens in 73 % of selected five-jet WW events. 

In all channels, the fitted W mass value is extracted using an unbinned maximum likeli- 
hood fit to the distribution of reconstructed 5C fit masses m in the region 70 < m < 88 GeV. 
The fit function is chosen empirically and consists of two terms: S{m) describes the signal 
contribution and B{m) the combinatorial and non-WW background. In the qqiv channel, 
the signal function consists of an asymmetric relativistic Breit-Wigner function with different 
widths above and below the peak: 

^^""^ = ^(m2-m2)2 + m2r2 2 ' ^^^^ 

where mo is the fitted mass and ri^2 are fixed parameters, Fi being taken for m < mo and r2 
otherwise, and A is a normalisation constant. In the qqqq channel the signal function S{m) is 
additionally multiplied by a Gaussian function of mean mo and width a, since this is found to 
improve the description of the reconstructed 5C fit mass distribution. These parameterisations 
were found to give adequate descriptions of the reconstructed distributions in Monte Carlo 
simulated data samples of around ten times the data luminosity. The parameters Fi and F2 
were determined using large samples of Monte Carlo signal events with mw = 80.33 GeV, 
and were parameterised as linear functions of ^/s. The parameter a for the qqqq channel was 
similarly determined and found to be independent of s/s. 

The contributions from combinatorial WW background, ZZ and Z/7 final states are rep- 
resented by the background function B{m), derived from Monte Carlo simulated events sep- 
arately at each centre-of-mass energy. The fractions of background assumed in the fits are 
fixed to those observed in Monte Carlo. As for the convolution fit, the fitted mass mo must be 
corrected for biases arising from initial-state radiation, the event selection, reconstruction and 
fitting procedures. Studies using Monte Carlo samples with the full range of ^/s values and 
true W masses from 79.33-81.33 GeV show these biases to have magnitudes of up to 500 MeV 
(in the qqrz/ and qqqq channels) , to be independent of -^/s, but to depend slightly on the true 
W mass. They were therefore parameterised as linear functions of the fitted mass using large 
Monte Carlo samples of both signal and background events, and applied as corrections to the 
raw fitted mass values. 

The Breit-Wigner fit is applied separately to the data taken at each centre-of-mass energy 
from 183 GeV to 209 GeV (dividing the 1999 and 2000 data samples into four and two energy 
bins respectively), and the results combined. The results for each channel (including both the 
modified P2.5 and Jo jet direction reconstruction methods in the qqqq channel) are shown in 
Table 2, together with the expected statistical errors evaluated using Monte Carlo subsample 
tests. The results are also shown as a function of data-taking year as the 'BW' points in 
Figure 4. Data taken in 1996 at ^/s « 172 GeV have not been reanalysed, and the Breit- 
Wigner fit results from [5] are shown. The reconstructed 5C mass distributions for qqiv 
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events and selected jet assignment combinations in qqqq events are shown in Figure 7(c) 
and (d), together with the fitted functions used to extract the W mass. 

8 Systematic uncertainties 

The main systematic uncertainties in the measurements of the W mass and width arise from 
the understanding of the detector cahbration and performance, the hadronisation of quarks 
into jets, possible final-state interactions in the qqqq channel, the modelling of non-WW 
background, the simulation of photon radiation in WW events and uncertainties in the LEP 
beam energy. These and other small systematic effects have been calculated separately for 
the qqiv and qqqq channels, using all three analysis techniques for the W mass, and for the 
convolution and reweighting fits for the W width. The determination of all systematic errors is 
described in detail below, and the results are summarised in Tables 4 and 5. Detector-related 
effects tend to increase slightly with energy; other uncertainties are taken to be constant 
unless stated otherwise. The magnitudes of the systematic uncertainties are generally rather 
similar between the different fitting techniques, but there are some significant differences, and 
these are also discussed below. 

8.1 Detector calibration and simulation 

The Monte Carlo descriptions of the jet and lepton energy scales, and energy and angular 
resolutions, are checked using samples of Z ^ qq, Z e+e~ and Z events recorded 

at ^/s « 91 GeV at the beginning of each data-taking year, and at other times during the 
data-taking periods in 1998, 1999 and 2000. Detailed comparisons between these data and 
corresponding Monte Carlo simulations are used to derive small adjustments at the level of 
reconstructed jets and leptons, which are then applied to the WW and background samples 
that are used to calibrate the bias corrections in the convolution and Breit-Wigner fits and to 
derive template distributions in the reweighting fit. Where necessary, the energy and angular 
resolution in the Monte Carlo simulation are degraded by applying Gaussian smearing to 
the reconstructed quantities. The associated systematic uncertainties are propagated to the 
measured W mass and width by varying the applied corrections. 

The properties of jets are checked using samples of two-jet and three-jet Z ^ qq events, 
and Z/7 — qq events from the high-energy LEP2 data samples, as follows: 

Jet energy scale: This is checked using Z ^ qq events reconstructed as two jets with the 
Durham algorithm and satisfying j/23 < 0.02. The same particle selection requirements 
and energy doTible-counting correction procedure are applied as for the WW analysis. 
The mean of the sum of the two jet energies is studied as a function of cos 9 = ^{\ cos 9i \ + 
I cos02|)) where 9i and 6*2 are the reconstructed polar angles of the two jets. The ratio 
of these energy sums in data and Monte Carlo is shown in Figure 8(a), and is used to 
derive corrections to the Monte Carlo energy scale as functions of jet cos 9 and data- 
taking year. The corrections in the forward region beyond cos 9 = 0.85 are much larger 
than in the central region, due to the difficulties in accurately modelling the complex 
detector geometry and larger amount of dead material. The residual uncertainty on 
the jet energy scale is 0.4%, dominated by contributions from Z data statistics, possible 
quark-fiavour dependences (assessed by repeating the studies after removing events with 
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Table 4: Summary of systematic uncertainties (in McV) on the measurements of the W mass. 
Results are given separately for the qq£z^ and qqqq channels (p2.5 jet direction reconstruction) 
with the convolution, reweighting and Breit-Wigner fitting methods. Results are also given 
for the convolution fit for the Jq and k_o.5 jet direction reconstruction methods, and for 
the combination of convolution qq£i/ and qqqq (P2.5) results (where the combination takes 
the systematic uncertainties and their correlations into account). The last column gives the 
systematic uncertainties for the mass difference Am(p2.5, /^-o.s) discussed in Section 9.1. 

a reconstructed secondary vertex indicating a heavy quark decay [28]) and possible 
variations during the course of a year. 

Jet energy resolution: The width of the distribution of two-jet energy sums is sensitive 
to the jet energy resolution, and was studied using the same techniques as the energy 
scale. The ratio of widths seen in data and Monte Carlo is shown in Figure 8(b) — the 
Monte Carlo resolution is about 4% too good for cosO < 0.85, and up to 20% too good 
in the forward region beyond cos^ = 0.85. After correction, the residual uncertainty 
lies between 0.6% and 2% depending on cos^, limited by Z data statistics. 

Jet energy linearity: The studies with two-jet Z — qq events check the energy scale for 
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Figure 8: Determination of energy corrections for jets (see text). Ratios of data to Monte 
Carlo are shown averaged over all data-taking years for: (a) jet energy scale as a function of 
cos 6*, (b) jet energy resolution as a function of cos9, (c) jet energy scale as a function of the 
jet energy itself. The results using the uncorrected simulation are shown by the filled points, 
and those with the corrected simulation are shown by the open points, with the error bars 
indicating the statistical error in each case. The horizontal lines indicate ratios of unity, and 
the dotted line in (c) shows the linearity correction used to parameterise the jet energy scale 
dependence on the jet energy itself. 
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Table 5: Summary of systematic uncertainties (in MeV) on the measurements of the W 
width. Results are given separately for the qq^z^ and qqqq channels with the convolution and 
reweighting fitting methods. The systematic uncertainties for the combination of <^tv and 
qqqq convolution fit results are also shown. 



^ 45 GeV jets, close to the average energy of jets produced in W decays, but event- 
by-event the latter range from about 20 GeV to 85 GeV. It is therefore important to 
check the linearity of the energy response, i.e. the energy scale for lower and higher 
energy jets. This has been studied by looking both at Z ^ qqg three-jet events and high 
energy Z/7 qq two-jet events. Coplanar three-jet events are selected by requiring 
y23 > 0.02 and 2/34 < 0.005, and that the sum of the inter-jet angles exceeds 355°. The 
jet energies can then be computed using the measured jet angles and masses, and the 
ratio of reconstructed to expected energies determined as a function of expected energy. 
The ratio of this quantity in data and Monte Carlo is shown in Figure 8(c), from which 
it can be seen that the energy scale in data is around 0.5 % higher for 30 GeV jets than 
for 45 GeV jets. 
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The behaviour at high jet energies is studied with Z/7 qq events taken at y/s > 
180 GeV, and satisfying 1/23 < 0.02 and \J s' / s > 0.85 where the reconstructed c+c~ 
colhsion energy after any initial-state radiation ^/s' is calculated as in [9]. In these 
events, the behaviour of the jet energy scale as a function of cos 6 is consistent with 
that seen for 45 GeV jets, but the overall energy scale is shifted downwards by about 
1 %, as can be seen for the high jet energy points in Figure 8(c). 

These studies are consistent with a linear dependence of the jet energy scale on the jet 
energy itself, with a slope of (—2.00 ± 0.30) x 10~^. The corresponding correction is 
applied to the energy scale in Monte Carlo simulation. The uncertainty is dominated by 
data statistics (0.26 x 10~^), but also includes systematic contributions from two-photon 
(0.06 X 10"^) and r-pair (0.09 x 10 ^) background modelling in the high energy qq sam- 
ples, and possible quark flavour dependences (0.09 x 10~^). Effects from hadronisation 
and four-fermion background modelling are found to be negligible. This uncertainty on 
the correction contributes a systematic error of around 4 MeV in the qqiu and 2 MeV 
in the qqqq W mass measurements. 

Although the data are consistent with a linear slope, a second order polynomial is 
also fitted and used to correct the simulation as an alternative. The corresponding W 
mass uncertainties when the curvature is varied within the range allowed by the data 
are 8 MeV and 2 MeV in the qqtu and qqqq channels. The final jet energy linearity 
uncertainties on the W mass and width are calculated as the quadrature sum of the 
shifts resulting from changing the linear correction by its uncertainty, and using the 
alternative second order polynomial correction with the maximum curvature allowed by 
the data. 

Jet angular resolution: The jet cos 9 and ^ resolutions are checked by using the two-jet 
Z — qq sample and studying the widths of the distributions of cos 0i-|-cos 62 and ^1 — ^2- 

These are found to be 4 % and 1 % narrower in Monte Carlo than data for the Jq jet 
direction reconstruction method, independent of cos^, and are smeared accordingly. 
The corresponding uncertainties are 0.4 % for cos 9 and 0.3 % for (f), dominated by Z data 
statistics. The uncertainties for the P2.5 direction reconstruction method are similar, 
though modelling of the jet angular resolution is somewhat worse. The differences 
between data and Monte Carlo are around a factor two larger than for the Jq direction 
reconstruction, necessitating correspondingly larger Monte Carlo corrections. 

Jet angular bias: A bias in the jet cos 6 reconstruction, equal in magnitude but opposite in 
sign for positive and negative cos 9, would not show up in the jet acollinearity measured 
by cos^i -I- cos ^2, and could have significant effects on the W mass and width. This is 
studied by using individual jets in two-jet Z ^ qq events, and calculating their cos 9 
separately using tracking and calorimeter information (the two detector systems have 
independent and uncorrelated angular reconstruction uncertainties). Some differences 
are seen, but these are generally well modelled by the Monte Carlo simulation. Since 
the jet cos 9 information is determined from both the tracking and calorimeter infor- 
mation, half the residual difference between data and Monte Carlo tracking-calorimeter 
deviations is taken as a systematic uncertainty on the absolute cos 9 bias. This study 
was performed for both the unmodified and alternative jet direction reconstruction al- 
gorithms, and no significant differences were seen. 
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Jet mass scale: No useful information on jet masses can be obtained from studies of Z 
data alone, and there are significant differences between the predictions of the Jetset, 
Ariadne and Herwig Monte Carlo models. The jet masses predicted by the Monte 
Carlo are therefore left unchanged by default, and a systematic error is assessed by 
scaling and smearing them event-by-event by the same factor as the corresponding jet 
energies — this is appropriate for the extreme case of jets composed entirely of massless 
particles. As an additional cross-check, the energies of all particles were scaled by the 
same factor as the jet energies, but before calculating the jet invariant masses. This gave 
results which were negligibly different from those obtained by scaling the jet masses. 

A similar approach is taken to study the modelling of leptons (electrons and muons sep- 
arately), using Z — e'^e" and Z — /x+/x~ events as follows: 

Lepton energy scale: The lepton energy scale is studied as a function of cos^ using the 
means of distributions of lepton energies. For electrons, the ratios of data to Monte 
Carlo means are typically within 1 % of unity, and are used to correct the Monte Carlo 
simulation, with a residual systematic uncertainty of 0.3 % including contributions from 
possible time dependence, the comparison of two independent event selections and data 
statistics. The ratios of data to Monte Carlo for muons typically agree to better than 
0.3 %, and a systematic uncertainty of 0.3 % is assigned, dominated by Z data statistics 
and the comparison of event selections. 

Lepton energy resolution: Studies of the width of the lepton energy distributions show 
that the electron energy resolution is around 19 % worse in data than Monte Carlo in the 
barrel region and 7 % worse in the endcap, and that the muon resolution is 6 % worse in 
data. These corrections are applied to the Monte Carlo with corresponding uncertainties 
of 2 %, dominated by data statistics. There arc also tails in the data resolution which 
are not well modelled by the Monte Carlo. For the W mass, a systematic error is 
conservatively estimated by doubling the Monte Carlo resolution correction, whilst for 
the width (which is much more sensitive to such tails) a more elaborate two-component 
smearing procedure is used to model both the core and tail resolution, with an additional 
uncertainty relating to the choice of smearing parameters. 

Lepton energy linearity: Possible dependences of the lepton energy scale on the lepton 
energy itself arc studied using e"'"e~7 and jji^ events taken both during the Z cali- 
bration and high energy running, by comparing the measured lepton energies with those 
determined from the track and cluster angles. No significant effects are seen, within a 
statistical precision on the slope of 3 x 10^^ for electrons and 6 x 10~^ for muons, 
and these values are used to assess the corresponding uncertainties on the W mass and 
width. 

Lepton angular resolution: The lepton cos 9 and (j) resolutions arc studied using the dis- 
tributions of cos^i -|- cos ^2 and (j)i — 02 in the same way as discussed above for jets. 
No evidence for Monte Carlo mis-modelling is seen, and the corresponding uncertainties 
are obtained from the statistical precision of the tests of 5-10 %. 

The effects of all these uncertainties on the various W mass and width analyses are shown 
in Tables 4 and 5. Uncertainties which affect only the qqez^ or qq/xz^ sub-channels are given 
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in terms of their effect on the combined qqfM results, and the uncertainties due to the mea- 
surements of jets and leptons are assumed to be uncorrelated. For the W mass, the most 
significant uncertainties are those associated with energy scales, jet angular biases and jet 
masses, whilst for the W width, energy resolution uncertainties play a bigger role. The uncer- 
tainties for the three fit methods axe rather similar, except for the jet energy scale uncertainty 
in the qqii' channel which is significantly smaller in the reweighting and Breit-Wigner fits 
than in the convolution fit. This is due to the different sensitivities of 4C and 5C kinematic 
fits with and without lepton information to variations of the jet energy scale, and the different 
use made by the three analysis methods of the different types of fit. 

A further correction is applied for the effects of beam-related background and detector 
noise, which lead to additional non-simulated occupancy in the detectors. This correction 
is evaluated by superimposing data events taken with a random beam-crossing trigger onto 
Monte Carlo WW and background events. The most significant effects come from additional 
clusters in the calorimeters (especially the hadron calorimeter in the forward region), and 
lead to shifts of around 10 MeV in the W mass and 2 MeV in the W width, with systematic 
uncertainties which are negligible in comparison to other detector-related effects. The effect 
of this noise on the data vs. Monte Carlo comparisons discussed above was also checked and 
found to be negligible. 

8.2 Hadronisation in W — > qq decays 

Uncertainties due to hadronisation in W — > qq decays are studied using large Monte Carlo 
samples where the same original four-fermion events have been hadronised using various 
different Monte Carlo models (string, colour dipole and cluster, as implemented in Jetset, 
Ariadne and Herwig respectively) and parameter sets (sec Section 2). These models have 
all been tuned to give a reasonable overall description of OPAL or ALEPH Z ^ qq data. 
Different models and tuned parameter sets describe particular features of the data to a greater 
or lesser extent, reflecting partly the emphasis placed on various variables {e.g. event shapes, 
charged and neutral particle multiplicities and fragmentation functions) by the different tune 
procedures. 

All these models give adequate descriptions of general event properties in W+W~ — > qqiv 
and W+W" qqqq events, and the limited data statistics do not allow any of the models to 
be disfavoured or excluded. However, they predict different fit biases or reweighting template 
distributions for the W mass and width fits. This is illustrated in Table 6, which shows the 
biases in fitted W mass and width from the convolution fit analysis applied to the same event 
samples hadronised with various different models, but calibrated using Monte Carlo events 
hadronised using JT. The statistical errors on the mass differences are calculated using a 
Monte Carlo STibsample technique and take into account the correlation between the samples 
due to the common initial four-fermion events. Taking the 'raw' mass shifts from Table 6, it 
is clear that the different models and tunes predict significantly different W mass biases, of 
up to ~ 40 MeV in both qqli/ and qqqq channels, and corresponding width biases of up to 
80 MeV. 

The mass biases have also been studied at the 'hadron' level, performing the jet finding 
on all stable hadrons^ produced by the Monte Carlo hadronisation model (before detector 
simulation), and then applying the full convolution analysis to these jets, together with the 

''Following the convention for Z decay multiplicities in [29], all particles with lifetimes of more than 3x 10"^" s 
were considered stable. 



34 



Models 




Raw mass 


shifts 






Adjusted mass shifts 






qqlv 


qqqq 


qqqq 


qqqq 
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qqqq 






P2.5 


Jo 


K-0.5 




P2.5 
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K-0.5 


JT'-JT 


-32 ±4 


-32 ± 4 - 


-40 ±3 


-33 ±5 


1±4 


4± 4 


5±3 


15 ±5 


AR-JT 


-25 ± 3 


-28 ±4 - 


-31 ±4 


-35 ± 5 


-7 ±4 


-6 ±4 


-4 ±4 


-6 ± 5 


AR'-JT 


-7 ±4 


-29 ±5 - 


-12 ±4 


6± 5 


1±4 


-20 ± 5 


2±4 


18 ± 5 


HW-JT 


-15 ±4 


6±4 


-3 ±3 


-15 ±5 


-13 ±4 


1±4 


-3 ±3 


-15 ±5 






Raw width shifts 












JT'-JT 


-35 ±9 




68 ±7 












AR-JT 


-43 ±9 




-9± 7 












AR'-JT 


-77 ±9 




-9 ±8 












HW-JT 


-5 ±9 




-42 ±8 













Table 6: Differences in W mass and width biases (in MeV) evaluated using the convolution fit 
for various different hadronisation models in comparison to the default JT model. For the W 
mass, results are given for the qq^i/ channel, and for the qqqq channel using the P2.5, Jo and 
K_0.5 jet direction reconstruction methods. Results are also given after adjusting the non-JT 
simulated events to have the same kaon and baryon content as that predicted by JT. The 
uncertainties are due to finite Monte Carlo statistics, and the definitions of the models JT, 
JT', AR, AR' and HW are given in Section 2. 

reconstructed leptons in the qq^i^ channel. The results show no significant differences between 
hadronisation models and tunes, showing that the biases are produced by the interplay of 
hadronisation and detector effects. At the hadron level, the various models predict different 
average jet masses for the two jets produced in W — >^ qq decays, but these differences are 
compensated by different inter-jet angles once the decays are boosted into the laboratory 
frame, where the invariant mass mjj of the jet-jet system is given by: 

ml = ml + ml + 2EiE2{l- /?i^2 cos ^12) (14) 

where nii, Ei and Pi are the mass, energy and velocity of jet i, and Ou is the angle between 
the two jets. The resulting average jet-jet invariant masses are therefore equal in all models. 
However, the dctcctor-level jet reconstruction introduces biases in both jet mass and angular 
distributions, and these biases are different in the various models, spoiling the hadron level 
compensation of the jet mass differences by the inter-jet angle differences and leading to 
significant differences in fitted W mass between models. 

A large part of the jet mass and angle biases is found to result from deficiencies in the 
reconstruction of kaons and baryons. In jet reconstruction, all charged particles are assigned 
the pion mass, and all neutral clusters zero mass. Charged kaons and protons (having m > 
m7r) will therefore be incorrectly reconstructed, as will and neutrons which in addition 
tend to have their energies badly estimated in the calorimeters. Although the hadronisation 
models have been tuned to Z data including kaon and baryon rates, there are significant 
differences between them, with e.g. JT and AR underestimating the production rates of 
kaons, and JT' and AR overestimating the production rates of baryons. Table 6 also shows 
the W mass differences between models after adjusting all the alternative models to have the 
same kaon and baryon content as JT.^ It can be seen that, particularly for the qqiv and 

^ Since particle multiplicities increase with the mass of the decaying boson, the event-by-event true W masses 
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qqqq channels with the Jq jet algorithm, the differences between models are greatly reduced 
after this adjustment procedure. Many other variables have also been studied, but no other 
significant dependences have been found, and the remaining mass differences are therefore 
taken to be indicative of genuine differences between the hadronisation models. 

In the qqqq analyses with modified jet direction reconstruction, significant mass bias dif- 
ferences persist after adjusting the kaon and baryon multiplicities. This is not surprising, as 
these reconstruction methods calculate the jet masses using one set of particles, and remove 
some particles or weight them differently when calculating the jet angles. This will again 
tend to spoil the cancellation between jet mass and angle differences seen at hadron level (see 
Equation 14), and the differences represent genuine uncertainty due to the modelling of hadro- 
nisation. The Monte Carlo models have been tuned to reproduce inclusive Z event properties, 
and it is not obvious that they can also be relied on to model e.g. event shape distributions 
when some particles are removed. To check this, event shape distributions were calculated 
in Z ^ qq events using only particles satisfying p > 2.5 GeV. The agreement between data 
and the various Monte Carlo hadronisation models for these distributions was found to be 
reasonable, and not significantly worse than for the inclusive distributions. Therefore, the 
models can also be expected to give a reasonable description of jet properties when particles 
with p < 2.5 GeV are removed. 

The final uncertainty on the W mass from the hadronisation of signal WW events is made 
up of two parts: the residual differences between hadronisation models after adjusting them 
to the same kaon and baryon production rates (where the largest difference between JT and 
any of the other models is taken), and an uncertainty related to the knowledge of kaon and 
baryon production rates in W — > qq decays. The latter are calculated from the measured K"*" , 
Kg and proton production rates in Z decays of 2.242 ± 0.063, 1.025 ± 0.013 and 1.048 ± 0.045 
[29]. The Kl rate is assumed to be equal to the Kg rate and the neutron to proton ratio to 
be that predicted by JT (0.97). The ratios of kaon (K"*" and K^) and baryon (proton and 
neutron) production in W and Z decays, R = n-w/nz, are taken to be Rk = 0.90 ± 0.01 and 
Rb = 0.95 ± 0.06, where the central values are taken from JT and the errors from the largest 
difference between JT and any other model.^ The measurements and ratios were combined 
to give predicted production rates of kaons and baryons in W decays of nK(W) = 2.94 it 0.06 
and ?7,b(W) = 2.00 ± 0.12. The uncertainties on these rates give small additional systematic 
errors on the W mass and width due to the dependence of the biases on kaon and baryon 
production. Small corrections are also applied to the fitted values to compensate for the 
differences between these predictions and those of JT. 

The differences in W width bias predicted by the various models are not related to kaon 
and baryon multiplicity, and do not reduce when the adjustment techniques used for the W 
mass are applied. No other variables have been found that play a similar role to the kaon and 
baryon multiplicities for the W mass. The largest raw difference between JT and any of the 
other models is therefore used to set the systematic uncertainty for the W width. 

Taking everything into account, the final uncertainties due to hadronisation in signal 
WW events are given in Tables 4 and 5. The systematic errors are dominated by the residual 
differences between hadronisation models, the uncertainties on kaon and baryon production 

and decay multiplicities are correlated, leading to 'artificial' W mass shifts when kaon and baryon multiplicities 
are adjusted. This efltect is removed using a second iteration in the adjustment procedure, which then changes 
the kaon and baryon content but leaves the average true W mass unchanged. 

^Kg, A and other hyperons are not included, since they typically decay into other particles which are already 
accounted for. 
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Mass shift (MeV) 




Width shift (MeV) 






P2.5 


P2.5 


P2.5 


Jo 




Jo 


Jo 


Model 




CV 


RW 


BW 


CV 


CV 


CV 


RW 


SK I {ki = 


0.9) 


19 ±2 


21 ±2 


16 ±2 


63 ±2 


113 ±3 


86 ±5 


80 ±3 


SK I {ki = 


2.3) 


41 ± 2 


41 ±2 


32 ±3 


125 ±2 


228 ± 3 


150 ± 5 


144 ± 5 


SK I (ki = 


100) 


136 ±3 


142 ±3 


105 ±4 


390 ±3 


674 ±4 


289 ±6 


295 ±6 


SK II 




-9 ±5 


24 ±6 


-5 ±3 


-3 ±5 


0±5 


33± 12 


26± 10 


SK II' 




-6 ±5 


-1±6 


-2 ±4 


-5 ±5 


4±5 


30± 12 


29± 10 


AR2-AR1 




29 ± 5 


27 ±6 


28 ±7 


66 ±5 


102 ±5 


128 ± 11 


104 ± 11 


AR3 




61 ± 25 


41 ± 24 


40 ±31 


145 ± 22 


251 ±31 


348 ± 49 


348 ± 25 


HW 




22 ±5 


35 ±8 


15 ±9 


42 ±5 


60 ±5 


27± 11 


55± 15 


LUBOEI BEC 


-24 ±8 


-22 ±9 


-27± 10 


-46 ±7 


-83 ±11 


41 ±14 


17 ±17 



Table 7: W mass and width shifts (in MeV) in the qqqq channel for various colour reconnection 
models and parameters, and for the LUBOEI Bose-Einstein correlation model, evaluated at 
the mean centre-of-mass energy of the data sample (196 GeV). Mass shifts are given for 
all three analysis methods with the p2.5 jet direction reconstruction method, and for the 
convolution fit with other jet reconstruction methods. The uncertainties are due to finite 
Monte Carlo statistics. 



contributing less than 5 MeV. The full analysis has been carried out for all three fitting meth- 
ods, which are found to have slightly different residual mass shifts and dependences on kaon 
and baryon multiplicities, particularly in the qq£i^ channel where the three methods employ 
4C and 5C kinematic fits in different ways. However, the overall hadronisation uncertainties 
are broadly similar. Uncertainties due to hadronisation in the Z/7 qq background are 
discussed in Section 8.5 below. 

8.3 Final-state interactions 

At LEP2 energies, the two W bosons produced in an e+e^ W+W~ event decay when their 
spatial separation is about 0.1 fm, much smaller than the typical hadronisation scale of 1 fm. 
The two hadronising W — qq systems in a W+W~ qqqq event therefore overlap, and their 
hadronisation may involve final-state interactions (FSI) between them, leading to exchange 
of four-momentum between the decay products of the two W bosons and possible biases in 
the reconstructed invariant mass spectra. Two sources of such interactions have been widely 
considered: colour reconnection and Bose-Einstein correlations (BEC). Both of these effects 
are well established in other systems, but are neither conclusively confirmed nor ruled out in 
W"'"W~ qqqq events [1]. The systematic uncertainties on the W mass and width from FSI 
effects are therefore determined by considering various phenomenological models which are 
consistent with the current limited knowledge of FSI in W^W" qqqq events. 

Colour reconnection effects in the perturbative phase have been shown to be small, giving 
rise to possible W mass biases of around 1 MeV [30]. However colour reconnection in the 
non-perturbative hadronisation phase may be substantial, and the mass and width biases can 
only be evaluated using Monte Carlo models. The models of Sjostrand and Khoze (SK I, 
SK II and SK II', as implemented in Pythia [30]), Lonnblad (AR2 and AR3 as implemented 
in Ariadne [31]), and Herwig [20] have been evaluated, and the resulting mass and width 
biases for all three analysis methods are summarised in Table 7. 
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The mass and width shifts for the SK I model arc evaluated using Monte Carlo samples 
generated with and without colour reconncction. The fraction of events with colour recon- 
nection depends on the SK I strength parameter kj. For each event, the reconnected or 
non-reconnected version is chosen according to a probability given by prec = 1 " exp(— FA;/), 
where V is the event-by-event space-time integrated product of the maximum colour field 
strengths of the two overlapping strings connecting the two quarks in each W ^ qq decay. 
The value of ki is not predicted by the model, and results for kj = 0.9, 2.3 and 100 are given 
in Table 7, corresponding to colour reconnection probabilities of 0.35, 0.57 and 0.97. Both 
the shifts and the reconnection probability for a given kj vary with centre-of-mass energy, 
e.g. from 37MeV at = 189 GeV to 47MeV at 207 GcV for kj = 2.3 in the convolution 
fit. They have therefore been evaluated at ^/s = 196 GeV, the mean centre-of-mass energy 
of the data sample, assuming linear dependences on -y/i. The convolution fit mass shifts for 
the P2.5 Jo and k-o.s jet direction reconstruction methods are also shown as functions of 
reconnection probability in Figure 9, together with the corresponding shifts for the W width 
measurement. The measured value of Am(p2.55 't-o.s) = —152 ± 68MeV (where the error is 
purely statistical — see Section 5.3) favours a reconnection probability of around 50 %; further 
discussion of this result is given in Section 9.1. Table 7 also gives results for the SK II and 
SK II' models, where the colour strings have infinitesimally small radii, but colour reconnec- 
tion occurs with unit probability on the first crossing (SK II) or only if it would reduce the 
total string length (SK II'). Both of these models predict smaller mass and width biases than 
SK I even with moderate values of kj, and are therefore not considered further when setting 
the systematic uncertainties on the W mass and width. 

Within the colour reconnection models implemented in Ariadne, reconnection occurs if 
it would reduce the total string length and is allowed within the constraints of the colour 
algebra factors [31]. In the first variant (ARl), colour reconnection occurs only amongst the 
decay products of one string (from a single W boson), and not between the two W bosons 
of a W"'"W~ — > qqqq event. In the AR2 model, colour reconnection is also allowed between 
strings, and hence between the two W bosons, but only for gluon energies below 2 GeV (the 
natural W width). The mass and width shifts due to colour reconnection are calculated as 
the difference between AR2 and ARl, to isolate the effect of colour reconnection between W 
bosons.^ Finally, in the AR3 model, colour reconnection is allowed between strings for all 
gluon energies, producing a strong effect and rather large mass and width shifts. However, 
this model is disfavoured both theoretically and by studies of three-jet events in Z data [33]. 
Additionally, studies of rapidity gaps in three-jet Z qqg events [34] disfavour the Ariadne 
implementation of colour reconnection even within one string. 

The Herwig Monte Carlo program also includes a colour reconnection model implemented 
in the framework of cluster hadronisation. In this model, a rearrangement of the association 
of partons to clusters occurs with a fixed probability of 1/9 if this rearrangement would lead 
to a smaller space-time extent of the clusters. The resulting mass and width shifts, shown in 
Table 7, are smaller than those for the SK I and ARiADNE-based models. 

All colour reconnection models studied show mass shifts which are reduced by factors of 
two to three by the P2.5 jet direction reconstruction method, and increased by factors of up 
to two for the K-0.5 method. Hence, as discussed in Section 4, the P2.5 method is used for 

*The 2 GeV limit on gluon energies for colour reconnection between W bosons is implemented in AR2 by 
running the dipole cascade twice, once down to 2 GeV with inter-W colour reconnection disabled, then again 
with it enabled. This results in some artificial additional high-energy showering, which was emulated in ARl 
for the purpose of this comparison by also running the cascade twice with an interruption at 2 GeV [32]. 
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Figure 9: (a) W mass and width shifts in the qqqq channel convolution fit as a function of re- 
connection probability prcc in the SK I model, for various different jet direction reconstruction 
methods (the points for the W width are slightly offset for clarity); (b) predicted mass differ- 
ence Am(p2.5, /t-o.s) ^ a function of reconnection probability Prec- The shifts are calculated 
at the displayed points using linear interpolation to y/s = 196 GeV, the mean centre-of-mass 
energy of the data sample, and the curves are drawn purely to guide the eye. The value and 
associated error (including both statistical and systematic contributions) measured from the 
data are indicated by the horizontal line and shaded band. 
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the main results of the qqqq channel analysis in this paper. The largest shifts are seen in the 
SK I model (depending on the value of fc/), and the AR2 model (AR3 being disfavoured both 
theoretically and experimentally). The final colour reconnection errors for the W mass and 
width are discussed in Section 9.1 below, where constraints from Am(p2.5, 1^-0.5) and studies 
of particle flow in W+W" qqqq events are also taken into account. 

Bose-Einstein correlations (BEC) between like-sign charged pion pairs are well established 
in both Z and W — > qq decays at LEP [35, 36, 37]. Although they are not implemented in the 
standard hadronisation models used in this paper, studies using dedicated samples generated 
using the LUBOEI [21] BEC model show that BEC between decay products originating from 
the same W boson (intra-W BEC) do not lead to significant W mass or width biases. However, 
just as in the case of colour reconnection, BEC between like-sign particles from different W 
bosons (inter-W BEC) may result in significant mass and width biases in the W+W~ — > qqqq 
analysis. These have been assessed using the LUBOEI BE32 model [21] as the difference 
between Monte Carlo samples generated with BEC affecting all possible like-sign particle 
pairs, and samples with BEC only between pairs from the same W boson. The parameters 
governing the properties of the generated correlations were tuned to describe BEC observed 
in Z — >^ qq decays as described in [36], and were set to A = 2.15 and R = 0.26 GeV. The 
resulting mass and width shifts are listed in Table 7. As in the case of colour reconnection, 
the mass shifts introduced by inter-W BEC are significantly reduced by the P2.5 alternative 
jet direction reconstruction method, and increased by the K-0.5 method. 

The possible existence of inter-W BEC has been experimentally investigated in [36, 37], 
and limits placed on its strength with respect to that predicted by the LUBOEI model. By 
fitting the BEC strength parameter A, the amount of inter-W BEC is measured in [36] to 
be a fraction 0.33 ± 0.44 of that predicted by LUBOEI, corresponding to a one standard 
deviation upper bound of 0.77. Assuming a linear relation between the BEC strength and 
the corresponding W mass and width shifts, the systematic errors on the W mass and width 
shown in Tables 4 and 5 are set to 77 % of those predicted by the LUBOEI model shown in 
Table 7. 

8.4 Photon radiation 

The dominant process contributing to four-fermion final states with an additional photon 

(e"'"e~ ffff7) is initial-state radiation (ISR) of photons from the incoming electrons and 
positrons, where the O(a^) treatment of KoralW is of more than adequate precision. How- 
ever, KoralW does not include all 0{a) photon radiation effects, e.g. radiation from the 
W bosons themselves (WSR), and interference between ISR, WSR and final-state radiation 
(FSR) from the outgoing charged leptons and quarks. A more complete treatment is provided 
by the so-called KandY [22] generator scheme, consisting of KoralW version 1.51 [22] and 
YFSWW3 [23] running concurrently. This introduces two major improvements over the Ko- 
ralW 1.41 samples used to calibrate the mass and width fits, namely the inclusion of 0{a) 
non-leading electroweak corrections and the screened Coulomb correction as opposed to the 
non-screened correction used previously. 

The differences in mass and width biases predicted by KandY and KoralW are shown 
for the convolution fit in Table 8. The results are similar for the other fit methods, and 
in the qqqq channel do not depend significantly on the choice of jet direction reconstruction 
method. As KandY gives a more complete treatment than KoralW, these shifts are applied 
as corrections to the final W mass and width results in this paper. However, the effects of 
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Mass (McV) 


Width (MeV) 




qqiv 


qqqq 


qqiv 


qqqq 


Shifts: 










KandY-KoralW 


-2 


1 


-22 


-22 


No non-leading EW correction 


17 


13 


14 


21 


No screened Coulomb correction 


-14 


-13 


15 


1 


ISR 0(a) - 0{a^) 


1 


1 


2 


1 


Uncertainties: 










Initial-state radiation 


1 


1 


2 


1 


Non-leading EW corrections 


8 


6 


7 


10 


Screened Coulomb correction 


7 


6 


8 


1 


Final state radiation 


1 


1 


2 


2 


Total uncertainty 


11 


9 


11 


10 



Table 8: Mass and width shifts measured using the convolution fit in the qqiu and qqqq 
channels for various changes to the treatment of photon radiation (see text) . The uncertainties 
due to finite Monte Carlo statistics are 1-2 MeV. The corresponding systematic uncertainties 
on the W mass and width are also given. 

the non-leading electroweak corrections and screened Coulomb corrections on the W mass 
partially cancel, and so are considered separately when assessing the total systematic error 
due to photon radiation. 

The main effect of non-leading electroweak corrections is to modify the ISR spectrum due 
to ISR-WSR interference [38], leading to bias in the W mass and width analyses, since ISR 
photons are not explicitly reconstructed in the kinematic fits. Studies in [38] show that the 
amount of ISR-WSR interference can be inferred from the rate of W"'"W~7 production with 
photons at large angles to the beam direction, and that the data arc described best by a 
parameter k = 0.38 it 0.47, where k = corresponds to the treatment in KandY and k = 1 
to that in KoralW. The data are therefore consistent with the prediction of KandY, and 
the uncertainty of 0.47 is used to determine the systematic uncertainties due to non-leading 
electroweak corrections. These are taken to be a fraction 0.47 of the mass and width shifts 
induced by reweighting KandY events to remove the corrections (see Table 8). 

The mass and width shifts induced by degrading the screened Coulomb correction of 
KandY to the unscreened correction implemented in KoralW are also shown in Table 8. 
The systematic uncertainties are taken to be half of these shifts. Finally, the systematic 
uncertainty due to the modelling of ISR is determined by reweighting events to degrade the 
0{a^) treatment of both KandY and KoralW to 0{a). The uncertainties due to FSR 
modelling have been assessed by reweighting events to change the rate of FSR from leptons 
by ±15%, and from quarks by ±50%, based on studies of Z and W decay data. The mean 
energies of FSR photons have also been varied by ±50%. The resulting uncertainties are 
very small, and also shown in Table 8. The total uncertainties due to photon radiation are 
determined from the quadrature sum of all the above sources, and amount to around 10 MeV 
for both the W mass and width. 
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8.5 Background 

In qqlu events, non-WW background is very small, except in the qqrz/ channel (see Table 1). 
The uncertainty due to modelling of background from Z /'j events is assessed by using Monte 
Carlo samples generated with P YTHIA (which has a simpler treatment of ISR) instead of KK2f, 
and by hadronising the KK2f samples with various different hadronisation models, as for signal 
WW events. The absolute rates of Z/7 and ZZ background are varied by ±20 %, following the 
uncertainties derived in [25]. The assigned uncertainty for ZZ background is larger than that 
for the on-shell Z-pair production cross-section [39], but includes contributions from other 
e"'"e~ff production diagrams. The errors assigned for all fit methods are given in Tables 4 
and 5. The total background error results from approximately equal contributions from 
hadronisation, Z/7 and ZZ background rate uncertainties. 

Background in the qqqq channel is much larger, and is dominated by e+e^ ^ Z/7 events 
giving a four-jet final state. Changing the hadronisation model used for the default KK2f Z/7 
samples from Jetset to Hekwig gives shifts of up to 20 MeV for the W mass and 32 MeV 
for the W width in the convolution fit, with Pythia and Ariadne lying in between Jetset 
and Herwig. The Monte Carlo modelling of this background is further investigated by using 
four-jet Z decays taken at ^/s « 91 GeV, and scaling the energies of all tracks and clusters 
by 200GeV/mz before applying the standard W+W~ qqqq event selection. These scaled 
events are then used in place of the standard Z/7 background samples in the W mass and 
width analysis. The mass and width shifts seen when using scaled Z events hadronised with 
Ariadne and Herwig instead of Jetset reproduce well the shifts seen in Z/7 events in all 
analysis methods, although there are some differences between the scaled Z and Z/7 mass 
distributions which are sensitive to details of the selection and scaling procedure. However, 
in all cases the scaled four-jet Z data events are found to lie between the predictions of 
Jetset and Herwig, so the differences between Jetset and Herwig Z/7 events at high 
energy are therefore used to set the systematic errors due to hadronisation in Z/7 events. 
Due to the small contribution of c+c~ ^ ZZ — >■ qqqq events to the selected sample (3-6 %) 
and the similar hadronic properties of W^W~ qqqq and ZZ — >■ qqqq events, the extra 
hadronisation uncertainty due to ZZ production is neglected. 

The absolute rate of Z/7 background events in the qqqq channel is varied by ±5 %, based 
on studies of the modelling of four-jet Z/7 background described in [25] and the modelling of 
the likelihood used in the W+W~ qqqq selection. The effect of changing the ISR modelling 
and ISR-FSR interference in KK2f is checked using the procedures described in [9] and found 
to be small. Finally, the absolute rate of ZZ — > qqqq events is varied by its uncertainty of 
±11% [39]. The total systematic errors due to the modelling of non-WW background are 
given in Tables 4 and 5, and are dominated by uncertainties in the modelling of fragmentation 
in Z/7 events. 

8.6 LEP beam energy 

Constraining the total reconstructed energy of the WW decay products to y/s greatly improves 
the event-by-event W mass resolution, but requires that the LEP beam energy be precisely 
known. The latter has been measured using a combination of resonant depolarisation and 
magnetic extrapolation based on NMR probes, complemented by measurements of the total 
bending field of the LEP dipoles using a flux loop, measurements of the beam energy using a 
dedicated spectrometer and studies of the accelerator synchrotron tune [40] . The beam energy 



42 



is known to a precision of between 10 MeV and 21 MeV, the largest uncertainty applying for 
the year 2000 data where special techniques were used to increase the beam energy to the 
highest possible value [40]. The beam energy uncertainties are largely correlated from year 
to year, and correspond to uncertainties of around 10 MeV on the W mass and 3 MeV on the 
width. 

Dispersion and other effects introduce a spread of between 160 MeV and 260 MeV in 
the event-by-event collision energy [40]. This effect, which is not included in the Monte 
Carlo simulations by default, introduces small shifts of around IMeV in the W mass and 
width, the full sizes of which are taken as additional systematic errors. A similar shift is 
produced by the average longitudinal boost of the collision centre-of-mass frame at OPAL of 
12-24 MeV [40]. This shift is caused by asymmetries in the distribution of the LEP radio- 
frequency accelerating system around the collider ring, and its full size is taken as an additional 
systematic uncertainty. Both these effects are included in the LEP beam energy entries of 
Tables 4 and 5. 

8.7 Other systematic errors 

Small discrepancies in the modelling of the data by the Monte Carlo lead to additional sys- 
tematic uncertainties in the convolution fit result. In both channels, the W mass and width 
bias depend weakly on the fitted event-by-event error, which on average is up to 1 % larger 
in data than Monte Carlo. Similarly, there are small discrepancies in the number of accepted 
jet assignment combinations in the qqqq channel (see Figure 1(c)). Both of these effects are 
assessed from the change in fit bias induced by reweighting Monte Carlo distributions to those 
of the data. There are no equivalent uncertainties in the reweighting and Breit-Wigner fits, 
where no significant discrepancies are seen in any important distributions. 

The bias corrections for the convolution and Breit-Wigner fits, and the template distribu- 
tions used in the reweighting fit, have small uncertainties due to finite Monte Carlo statistics. 
Such uncertainties also play a role in e.g. the hadronisation and detector systematics, but 
are accounted for in the corresponding systematic errors where appropriate. 

8.8 Consistency checks 

The complete W mass analysis is performed for the convolution, reweighting and Breit-Wigner 
fits, and the width analysis for the convolution and reweighting fits, with an additional 5C 
convolution fit in the qqiu channel. The results are given in Tables 2 and 3. The differences 
between the various fit results are summarised in Table 9, together with the expected RMS 
differences evaluated using a large number of common Monte Carlo simulation subsamples 
analysed by each method in a consistent manner. The observed differences between fit meth- 
ods are compatible with expectation. The consistency between the different fit methods, 
data-taking years and analysis channels can also be seen in Figures 4 and 5. Given these 
results, no additional systematic uncertainty associated with any individual fit method is as- 
signed. Complete systematic error analyses are also performed for all fitting methods except 
the CVS qqiu width fit (see Tables 4 and 5), and largely comparable uncertainties obtained, 
again giving confidence in the results. 

The statistical correlations between the fit methods have also been assessed using similar 
Monte Carlo subsamples, and found to be between 0.65 and 0.88. Combining the results from 
all three methods would only reduce the statistical error on mw by around 2 %. Given this. 
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Channel 




Fitted m-w 






Fitted F-yv 






CV-RW 


RW-BW 


BW-CV 


CV-RW 


RW-CV5 


CV5-CV 


qqei^ 


19 ±48 


-8 ±55 


-11 ±59 


-313 ± 170 


34 ± 160 


279 ± 160 




-56 ± 52 


-35 ± 56 


91 ±68 


35 ± 180 


8 ±160 


-43 ± 180 


qqTv 


65 ±92 


20 ±89 


-85 ± 116 


-326 ± 260 


-115 ±270 


441 ± 290 


qqii' 


0±33 


-6 ±35 


6 ±42 


-162 ± 100 


-15 ±90 


177 ± 90 


qqqq(p2.5) 


45 ±45 


22 ±56 


-67 ± 54 








qqqq(Jo) 


11 ±39 


-41 ± 43 


30 ±40 


-51 ±97 







Table 9: Differences in fitted W mass and width values between pairs of fit methods in each 
analysis channel. The uncertainties indicate the expected RMS spread in the fitted differences, 
evaluated using Monte Carlo simulation subsamples. 

and the large correlation between the systematic uncertainties of the different methods, they 
are not combined and the final results are taken from the convolution fit alone. 

9 Results 

The final results of the W mass and width fits are presented in this section, taking into account 
both statistical and systematic errors. The results for the measurement of Am(p2.5, /^-o.s) 
are given in Section 9.1 and used in conjunction with previous measurements to derive a limit 
on the SK I model parameter kj. The results for the W mass and width incorporating this 
limit are then given in Section 9.2. 

9.1 Colour reconnection limit 

The differences in W mass Am{ Jx, 1^-0.5) extracted using an alternative jet direction recon- 
struction method X and method K-0.5 are sensitive to possible final-state interactions in the 
W+W~ qqqq channel and can be used to set a limit on colour reconnection, as discussed 
in Section 5.3. Of all the direction reconstruction methods considered, Monte Carlo studies 
show that the mass difference Am(p2.5, /^-o.s) is the most sensitive to the SK I model, for 
both moderate and large amounts of colour reconnection, and is therefore used to set a limit 
on kj. 

The systematic uncertainties on Am(p2.5; '^-o.s) are evaluated using the same techniques 
as discussed in Section 8 and are given in Table 4. Detector effects are largely correlated 
between jet algorithms and the resulting residual uncertainties are very small. Hadronisation 
effects are calculated by taking the effect on Am(p2.5, /^-o.s) when simultaneously changing 
the hadronisation models used to set the bias corrections in both mass measurements, and are 
dominated by the difference between JT and AR' (see Table 6). Background and modelling 
uncertainties are assessed similarly. The largest uncertainty comes from Bose-Einstein cor- 
relations, which have a similar effect to colour reconnection on the value of Am(p2.5, 1^-0.5), 
albeit with the opposite sign. Taking all uncertainties into account, the result is: 

Am(p2.5, K-0.5) = -152 ± 68 ± 41 ± 45 MeV, 

where the first error is statistical, the second systematic excluding Bose-Einstein correlations 
and the third from Bose-Einstein correlations. Within the SK I model, this corresponds to a 
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ki 


Prec 


Uncertainty 






Mass 


Width 


0.9 


0.35 


19 


86 


1.5 


0.47 


31 


119 


1.8 


0.52 


35 


133 


2.3 


0.58 


41 


151 


3.1 


0.65 


50 


172 



Table 10: Systematic uncertainties (in MeV) on the qqqq channel measurements of the W 
mass and width in the SK I model for various values of the model parameter kj and associated 
colour reconnection probability prec- The uncertainties corresponding to kj = 2.3 are used for 
the final result. 

value of kj = 1.7l^'2, ^ colour reconnection probability of prec = 0.5i;i;o 27 \/s=196GeV 
(see Figure 9). 

This result is combined with the OPAL limit determined from studies of particle flow in 
inter-jet regions of W''~W~ qqqq events [41].^ Statistical correlations between the two 
methods were evaluated using Monte Carlo subsample techniques, and found to be less than 
10 %. The largest systematic errors from hadronisation and Bose-Einstein correlation were 
taken to be fully correlated. The final result is prec = 0-43lo20' with the one standard 
deviation upper bound of p^ec < 0.58 corresponding to ki = 2.3. This value of ki has 
been used to calculate the systematic errors due to colour reconnection for the W mass and 
width shown in Tables 4 and 5, based on the SK I model. The systematic errors which 
would result for other values of kj can be seen in Figure 9 and are listed in Table 10. The 
corresponding limit at 95% confidence level is Prec < 0.70 or kj < 4.0. The particle flow 
and Am(p2.5) '^-o.s) results are not sensitive enough to test the predictions of the AR2 and 
Herwig colour reconnection models, but these give mass and width shifts smaller than those 
predicted by the SK I model with kj = 2.3, as can be seen from Table 7. 

The final W mass and width fit results are not adjusted to compensate for the possible 
effects of colour reconnection. A constant colour reconnection uncertainty is assumed, inde- 
pendent of ^/s and evaluated at the mean centre-of-mass energy of the data sample. This 
avoids changing the relative weights of the different energy points as a function of the energy 
dependence of the colour reconnection uncertainty of any particular model. 

9.2 Results for the W mass and width 

The results for the W mass and width, including both statistical and systematic errors, are 
given for all analysis methods in Table 11. The central values differ slightly from those 
in Tables 2 and 3 as systematic errors have been taken into account in the combination of 
results from different years. Combining both qqiu and qqqq channels using a minimisation 
technique, and assuming all systematic errors to be correlated between the two channels, the 
final results (taken from the convolution fit and using the running width scheme for the 
Breit-Wigner distribution as implemented in KoralW) are: 

mw(qq^i^ + qqqq) = 80.416 ± 0.042 ± 0.031 ± 0.009 GeV , 

^The analysis of W decay charged multiplicity diflFerences between W'''W~ qqiv and W^W~ — > qqqq 
events presented in [41] is not sensitive enough to colour reconnection to make any significant difference to the 
combination presented here. 
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Convolution 


W mass (GeV) 
Reweighting 


Breit-Wigner 


qqlv 

qqqq(p2.5) 

qqqq(Jo) 


80.449 ± 0.056 ± 0.028 
80.353 ± 0.060 ± 0.058 

80. 394 ± 0.05^ ± 0.755 


80.451 ± 0.058 ± 0.022 
80.308 ± 0.064 ± 0.056 

80.383 ± 0.056 ± 0.136 


80.457 ± 0.063 ± 0.029 
80.278 ± 0.072 ± 0.057 

8O.4I6 ± ± C.i57 


Combined 


80.416 lb 0.042 ± 0.032 
Convolution 


80.405 ± 0.044 ± 0.028 

W width (GeV) 
Reweighting 


80.390 ± 0.048 ± 0.032 


qqiv 
qqqq(Jo) 


1.927 ±0.135 ±0.091 
2.125 ±0.112 ±0.177 


2.088 ±0.131 ±0.085 
2.176 ±0.130 ±0.180 




Combined 


1.996 ± 0.096 ± 0.102 


2.113 ±0.101 ±0.097 





Table 11: Summary of W mass and width results for all fit methods. In each case, the 
first error is statistical and the second systematic, including the error on the LEP beam 
energy. The qqqq mass results with the Jq jet direction reconstruction algorithm are shown 

for comparison purposes, and are not included in the combination. The results quoted for the 
Breit-Wigner fit include the 172 GeV results from [5], as these data have not been reanalysed 
using this fit method. 

rw(qq^z^ + qqqq) = 1.996 ± 0.096 ± 0.102 ± 0.003 GeV , 

where in each case the first error is statistical, the second systematic and the third due to 
the uncertainty in the LEP beam energy. The estimated correlation between the two results 
is 0.04. The full breakdown of systematic errors is given in Tables 4 and 5. The qqqq 
channel has a weight of 0.34 in the mass measurement and 0.35 in the width measurement; 
the qqqq weights for the other fitting techniques are similar. In the absence of systematic 
uncertainties due to final-state interactions, and using the Jq jet direction reconstruction 
method, the statistical error of the combined W mass measurement would be 0.038 GeV, only 
10 % smaller than the 0.042 GeV of the present result. This demonstrates that the modified 
jet direction reconstruction technique significantly reduces uncertainties due to final-state 
interactions, and allows most of the statistical power of the qqqq channel to be exploited. 
The differences between the fitted values of mw and Fw in the qqqq and qq£u channels 

are: 

Amw (qqqq - qqM = -0.097 ± 0.082 ± 0.039 GeV , 
Arw(qqqq - qqM = 0.198 ± 0.175 ± 0.124 GeV , 

where in each case the first error is statistical and the second systematic, excluding uncer- 
tainties due to possible final-state interactions in the qqqq channel. For these results, the 
hadronisation and background uncertainties are conservatively taken to be uncorrelated be- 
tween the qq^z^ and qqqq channels, and all other uncertainties are taken to be fully correlated. 
Significant non-zero values of Amw or AFw could indicate that final-state interactions are 
biasing the values determined from the qqqq channel; however these values are consistent 
both with zero and with the shifts expected from colour reconnection in the SK I model with 
ki = 2.3. 

The result for the W mass is combined with previous OPAL measurements using W"*" W~ — >■ 
ii/lu events at ^/s values between 183 GeV and 209 GeV [3] and from the dependence of the 
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WW production cross section on mw at -y/i 161 GeV [4]. The final OPAL result for the W 
mass is then 

mw = 80.415 ± 0.042 ± 0.030 ± 0.009 GeV , 

where again tlic first error is statistical, the second systematic and the third due to uncer- 
tainties in the LEP beam energy. 

10 Conclusions 

The mass and width of the W boson are measured using e"'"e~ W^W~ events from the 
complete data sample collected by OPAL at centre-of-mass energies between 170 GeV and 
209 GeV, using event-by-event reconstruction of the W mass in the qq£u and qqqq final 
states. The result for mw is combined with earlier OPAL results using iuiv events [3] and 
the dependence of the WW production cross-section on mw at threshold [4]. The final results 
are: 

mw = 80.415 ± 0.042 ± 0.030 ± 0.009 GeV , 
Pw = 1.996 ± 0.096 ± 0.102 ± 0.003 GeV , 

where the first error is statistical, the second systematic and the third due to the uncertainty 
on the LEP beam energy. These results are consistent with, and supersede, our previous 
results [5, 6, 7], and are also consistent with other values measured at LEP [42] and the 
Tevatron [43]. The results are also consistent with the values inferred indirectly from pre- 
cision electroweak data [1] , providing a powerful consistency test of the Standard Model of 
electrowcak interactions. 

Limits arc placed on the strength of possible final-state interactions in e^e~ W^W~ — > 
qqqq events, by studying the evolution of the fitted W mass measured with various jet algo- 
rithms having differing sensitivities to such interactions. In combination with OPAL results 
based on particle flow in the regions between jets in c+c" W"'"W~ qqqq events, a 95 % 
confidence level upper limit on the SK I colour reconnection model parameter kj is set at 
kj < 4.0, corresponding to a colour reconnection probability of 0.70. 
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